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High-Tc Superconductivity and Antiferromagnetism in Multilayered Copper Oxides 
— A New Paradigm of Superconducting Mechanism — 

Hidekazu Mukuda*, Sunao Shimizu^, Akira lyo^-'-, and Yoshio Kitaoka^ 

Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531 
^National Institute of Advanced Industrial Science and Technology (AIST), Umezono, Tsukuba 305-8568 

High-temperature superconductivity (HTSC) in copper oxides emerges on a layered Cu02 
plane when an antiferromagnetic Mott insulator is doped with mobile hole carriers. We review 
extensive studies of multilayered copper oxides by site-selective nuclear magnetic resonance 
(NMR), which have uncovered the intrinsic phase diagram of antiferromagnetism (AFM) and 
HTSC for a disorder-free Cu02 plane with hole carriers. We present our experimental findings 
such as the existence of the AFM metallic state in doped Mott insulators, the uniformly mixed 
phase of AFM and HTSC, and the emergence of d-wave SC with a maximum Tc just outside 
a critical carrier density, at which the AFM moment on a Cu02 plane disappears. These 
results can be accounted for by the Mott physics based on the t-J model. The superexchange 
interaction Jin among spins plays a vital role as a glue for Cooper pairs or mobile spin-singlet 
pairs, in contrast to the phonon-mediated attractive interaction among electrons established 
in the Bardeen-Cooper-Schrieffer (BCS) theory. We remark that the attractive interaction for 
raising the Tc of HTSC up to temperatures as high as 160 K is the large Jin (~ 0.12 eV), which 
binds electrons of opposite spins to be on neighboring sites, and that there are no bosonic glues. 
It is the Coulomb repulsive interaction i7 (> 6 eV) among Cu-3d electrons that plays a central 
role in the physics behind high-Tc phenomena. A new paradigm of the SC mechanism opens to 
strongly correlated electron matter. 

KEYWORDS: copper oxide, correlated electrons, high-Tc superconductivity, antiferromagnetism, phase 
diagram, NMR 



1. Introduction 

In 1986, Bednorz and Miiller unveiled a new class of 
ceramic superconducting materials composed of a lay- 
ered structure of a two-dimensional Cu02 plane. ^ Subse- 
quently, the superconducting transition temperature (Tc) 
in cuprates rises up to '^133 K in HgBa2Ca2Cu308+5 
(Hgl223),2 which reaches rc=164 K under pressure. 
The discovery of a remarkably high Tc in the cuprate 
family has not only opened many possibilities for poten- 
tial technical applications, but has also provided a chal- 
lenging research subject for condensed-matter physics 
and material sciences. 

Despite more than a quarter of a century of research, 
there is still no universally accepted theory for the mech- 
anism of high-Tc superconductivity (HTSC) in cuprates. 
The main controversy arises with regard to the origin of 
the attractive force for the formation of Cooper pairs, 
which leads to such a remarkably high SC transition. In 
conventional superconductors with a relatively low Tc, 
the phonon-mediated electron-electron interaction is the 
attractive force for the formation of Cooper pairs in the 
Bardeen-Cooper-Schrieffer (BCS) theory, which was es- 
tablished half a century ago.^ The HTSC in cuprates 
emerges on a layered structure of a two dimensional 
Cu02 lattice when an antiferromagnetic Mott insula- 
tor is doped with hole carriers. The strong hybridiza- 
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tion of Cn-'idr^2 _yi and 0-2pa orbitals brings a large su- 
perexchange interaction Jjn ~ 0.12 eV (~1300 K) in the 
Cu02 plane for nondoped cuprates. ''^^'^ Therefore, an in- 
timate relationship between antiferromagnetism (AFM) 
and HTSC is believed to be a key for understanding the 
origin of the remarkably high SC transition in cuprate 
superconductors. "'^^"^'^ Experimentally, however, in pro- 
totype high-Tc cuprates La2-xSr^Cu04 (LSCO), AFM 
and HTSC are separated by the spin-glass phase. In 
LSCO, the chemical substitution of Sr^+ for La'^^ is nec- 
essary to increase the planar Cu02 hole density, but it 
introduces some local disorder simultaneously into nearly 
two-dimensional Cu02 planes and buckling on a CuOe 
octahedral unit, which make doped hole carriers local- 
ize owing to the Anderson localization mechanism. As a 
result, it is inevitable that the intrinsic electronic char- 
acteristics of LSCO are masked so that AFM and HTSC 
are separated by the spin glass phase. 

Multilayered cuprates provide us with the opportu- 
nity to research on the characteristics of a disorder- 
free Cu02 plane with hole carriers homogeneously 
doped. Figures 1(a) and 1(b) respectively show the 
crystal structures of n-layered cuprates in the se- 
ries of A/Ba2Ca„_iCu„02„+2+(5 (M=Hg, Tl, and Cu) 
and Ba2Ca„_iCu„02,i(Fj^Oi_j^)2, denoted as M12{n- 
\)n and 02(n-l)nF. Here, n is the number of Cu02 
planes within a unit cell. Copper oxides with more than 
three layers comprise inequivalent types of Cu02 lay- 
ers: an outer Cu02 plane (OP) in a five-fold pyrami- 
dal coordination and an inner Cu02 plane (IP) in a 
four-fold square coordination. Site-selective nuclear mag- 
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Fig. 1. (Color online) Crystal structure of n-layered cuprates of (a) A/12(ra-l)n and (b) 02(n-l)nF. Copper oxides with more than three 
layers comprise inequivalent types of Cu02 layers: an outer plane (OP) in a five-fold pyramidal coordination and an inner plane (IP) 
in a four-fold square coordination. Although the disorder may be introduced along with the chemical substitution in charge reservoir 
layers (CRLs), it is effectively shielded out of OPs, and hence homogeneously hole-doped Cu02 planes with ideal flatness are realized, 
especially at IP, which is ensured by the narrow NMR linewidths (see Fig. 4). 



netic resonance (NMR) and nuclear quadrupole reso- 
nance (NQR) studies are unique tools for differentiat- 
ing layer-dependent electronic characteristics microscop- 
ically. '^^"^'^ One of the remarkable features and advan- 
tages of multilayered cuprates is that the Cu02 layers 
are very flat and homogeneously doped, which have been 
ensured by the narrowest NMR linewidth to date among 
the very high quality cuprates investigated thus far (for 
example, see Fig. 4). In multilayered cuprates, the car- 
rier densities are inequivalent between OP and IP ow- 
ing to an imbalance in the Madelung potential at each 
Cu02 plane. Namely, since IP is farther from the charge 
reservoir layer (CRL) than OP, the carrier density at IPs 
is always lower than that at OPs. Carrier density can 
be tuned by the oxygen deficiency in CRLs {MOs) for 
M12(n-l)n or by the chemical substitution of F at apical 
oxygen sites for 02(n-l)nF. Such chemical substitutions 
introduce some disorder in CRLs, which may be slightly 
mapped onto OP, but the disorder potential at IP can 
be effectively shielded owing to the presence of conduct- 
ing OP, as deduced from the narrower NMR linewidth at 
IP than at OPs (for example, see Figs. 3 and 4). In this 
context, ideally flat Cu02 planes are realized especially 
at underdoped IP, differentiating multilayered cuprates 
from monolayered cuprate LSCO. 

In this paper, we review a decade of extensive NMR 
investigations of n-layered cuprates with n=3, 4, and 
5, which have revealed the intimate relationship be- 
tween AFM and HTSC for a disorder-free Cu02 plane 
with hole carriers homogeneously doped. The in- 
trinsic phase diagram possesses the following features: 
The AFM metallic state is robust and coexists uniformly 
with the HTSC at a single Cu02 plane in a region ex- 
tending up to the optimally doped one. The critical car- 
rier density pc at which the AFM order collapses de- 
creases from 0.10 to 0.08 to 0.075 as the interlayer mag- 
netic coupling becomes weaker when decreasing from 
n=5 to 4 to 3, respectively. This provides a reason- 
able explanation why the AFM order in n=l:LSC0 and 



n=2:YBa2Cu306-i-a;(YBC06-i-a:) collapses at carrier den- 
sities with pc=0.02 and 0.055, respectively. We reveal 
that the SC gap and Tc exhibit a maximum irrespec- 
tive of n at p ~ 0.16 just outside Pc{Mafm=0) ~ 0.14, 
where the AFM moment (Mafm) inherent in the Cu02 
plane totally disappears in the ground state. We high- 
light that the ground-state phase diagram of AFM and 
HTSC (see Fig. 22) is in good agreement with the 
ground-state phase diagrams in terms of either the t- 
J niodel,"-^^' 13-22, 24-28, 31, 32 ^^le Hubbard model in 

the strong-correlation regime. ^^'^^ The results presented 
here demonstrate that the in-plane superexchange inter- 
action Jin plays a vital role as a glue for Cooper pairs or 
mobile spin-singlet pairs, which will lead us to a coherent 
understanding why Tc is so high for hole-doped cuprates. 

2. Experimental 

2.1 Sample preparation and characterization 

Polycrystalline powder samples of n-layered cuprates, 
i.e., A/Ba2Ca^_iCuj^02n-t-2-t-5 

(M=Hg, Tl, and Cu) and 
apical F-substituted Ba2Ca„_iCu„02n(FyOi_y)2, were 
prepared by the high-pressure synthesis technique, as 
described in the literature. ^^"^^ To obtain more under- 
doped samples of the Hgl245 system, as-prepared sam- 
ples in nearly optimally doped region [Hgl245(OPT)] 
were annealed in a quartz tube with Cu powder or in 
Ar gas atmosphere for more than several hundreds of 
hours.39'46,62 Ba2Ca„_iCu„02„(Fj,Oi_j,)2, the sub- 
stitution of oxygen (O"^) for apical fluorine (F^^), i.e., 
a decrease in nominal fluorine content (y), results in the 
doping of holes into Cu02 layers, increasing Tc-^"^^^^ Al- 
though the real apical fluorine F"^ content may devi- 
ate slightly from the nominal one, Tc and carrier den- 
sity (p) can be tuned by changing the nominal content 
y in this series, which provides an opportunity to in- 
vestigate the characteristics of Cu02 layers over a wide 
p range systematically in the homologous series of n- 
layered cuprates. 

Powder X-ray diffraction measurements indicate that 
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Fig. 2. (Color online) (a) Lattice parameters a and c plotted 
against n for the Hgl2(n-l)n and 02{n-l)nF systems. Note 
that the c/2 of 02(n-l)nF is compared with the c of Hgl2(n-l)?i, 
because of the difference in the unit cell, (b) Relationship be- 
tween Tc and n for the homologous series of n-layered cupratcs 
M12(ra-l)n and 02(n-l)nF. [cited from refs.S'^'fOj 



the samples used for NMR/NQR measurements are al- 
most entirely composed of a single phase. As shown in 
Fig. 2(a), the c-axis length monotonically increases with 
increasing n, which can be fitted using linear functions, 
c{n) ~ 9.451 + 3.171 x (n - 1)[A] for Hgl2(n-l)n and 
c(n)/2 ~ 7.205 + 3.191 x (n - 1)[A] for 02(n-l)nF."'6O 
The first term corresponds to the distance between OPs 
through CRL averaged in n-layered cuprates; the dis- 
tances are 9.451 [A] for Hgl2(n-l)n and 7.205[A] for 02(n- 
l)nF, originating from the difference in the structure of 
CRL between the two systems. The coefficient of the sec- 
ond term corresponds to the average distance between 
adjacent Cu02 planes, which is almost equal to the c-axis 
length of the infinite layer CaCu02(IL) (c(oo)= 3.179 A). 
The a-axis length of two systems also approaches that 
of the CaCu02(IL) (a(oo)= 3.856 A) as n mcreases, as 
shown in the lower panel of Fig. 2(a). 

The values of Tc of all the samples were uniquely de- 
termined by susceptibility measurement using a SQUID 
magnetometer, which exhibits a marked decrease due to 
the onset of SC diamagnetism. The Tc of nearly opti- 
mally doped samples exhibits a striking dependence on 
n, as shown in Fig. 2(b), which has been confirmed in a 
homologous series.^^'^"'^^'^^ 

2.2 NMR/NQR measurements 

For NMR/NQR measurements, the powder samples 
were aligned along the c-axis at an external field (B) of 
~16 T and fixed using sty cast 1266 epoxy. Figure 3 shows 
a typical ^^Cu-NMR spectrum of n=5:Hgl245(OPT)tJl, 
in which the spectra from OP and IP are sepa- 
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^^Cu-NMR spectrum of 
®^Cu-NMR spectra at a pyramid- 
type outer Cu02 plane (OP) and a square-type inner one (IP) 
are separately observed owing to the differences in Knight 
shift. ^® The linewidths in the ^^Cu-NMR spectra are as narrow 
as 50 Oe at IP and 110 Oe at OP at ~15 T(_B||c), indicating 
that IPs are ideally flat and homogeneously holc-dopcd. [cited 
from ref.'^^l 



rately observed owing to the difference in Knight shift, 
which enables us to study multilayered compounds site- 
selectively. '^'^"^'^'^^"^^ Owing to the difference in local 
structure between IP and OP, the ^"^Cu-NQR frequen- 
cies {^^lyg) at OP and IP are typically '^16 and ~8.4 
MHz, respectively, for n=5:Hgl245(OPT)ttl,3S in which 
the presence of apical oxygen on OP generally makes 
the ^^i^Q of OP large. Moreover, it is remarkable that 
the linewidths in the ^•^Cu-NMR spectra are particularly 
narrow in multilayered cuprates, which are as narrow as 
50 Oe for IP and 110 Oe for OP in ?i=5:Hgl245(OPT)ttl 
even at ~15 T(B||c). For comparison, the ^^Cu-NMR 
spectra of typical n-layered cuprates are presented in Fig. 
4: (a) n=l : LSCO (x=0.24, Tc=18 K)^^ and Hgl201 
(Tc=96 K),69 (b) n=2: YBCO7 (Tc=92 K) and Bi2212 
(Tc=86 K),'^o (c) n=3: Hgl223 (Tc=133 K)^^ and Cul223 
(Tc=71 K),37 (d) n=4: Cul234 (Tc=:117 K),36.37 and (e) 
n=5: Hgl245(OPT)Ul (Tc=108 K).^* Here, the spectra 
are displayed against {B — B^cs) / Bj^os for the normal- 
ization of differences in field conditions, where the Bics 
is the resonance field for ^■^Cu. The broadening of ^'^Cu- 
NMR linewidth originates from the inhomogeneity of the 
Knight shift, and the distribution of the quadrupole shift 
when Cu02 planes are buckled. Thus, the figure indicates 
that the ^^Cu-NMR hnewidth becomes narrower as n 
increases, suggesting that multilayered cuprates possess 
very flat Cu02 layers from a microscopic point of view 
as well as homogeneous electronic states over the Cu02 
plane with hole carriers. Moreover, we also note that the 
SC in Hgl201(n=l) and YBC07(n=2) compounds also 
occurs on the Cu02 plane with less disorder, which is 
one of the key factors for the relatively high SC transi- 
tion at Tc=96 and 92 K among n=l and 2 compounds, 
respectively.''^ These facts that the Cu02 planes are ide- 
ally flat and homogeneously hole-doped in multilayered 
cuprates enable us to investigate the intrinsic properties 
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Fig. 4. (Color online) ^Scu-NMR spectra of (a) n=l: LSCO 
(x=0.24, Tc=18 K)68 a^nj Hgl201 (Tc=96 K),^^ (b) n=2: 
YBCO7 {Tc=92 K) and Bi2212 (Tc=86 K),™ (c) n=3: Hgl223 
(Tc=133 K)'^'^ and Cul223 (Tc=71 K),^'^' (d) n=4: Cul234 
(rc=117 K),36.37 a^nd (g) „^5. Hgl245(OPT)ttl (Tc=108 K).38 
Here, the spectra are displayed against {B — i?rcs)/Brcs for the 
normalization of differences in field conditions, where -Bros is a 
resonance field for ^'^Cu. These spectra were measured under the 
following experimental conditions; LSCO (-Bros -L c=10.79 T, 
T=20 K, aligned polycrystal{APC)),^** Hgl201 (Bres || c=7.46 
T, T=100 K, APC),'59 YBCO7 (Bres II c=15.08 T, T=90 K, 
single crystal), Bi2212 (Sres || c=10.95 T, T=77 K, single crys- 
tal),''" Hgl223 (Bres II c=10.95 T, r=140 K, APC),'''' Cul223 
(Bres II c=15.24 T, T=160 K, APC),^'' Cul234 (Bros II c=15.25 
T, r=90 K, APC),^'^'^^ and Hgl245 (Bros |j c=15.25 T, T=280 
K, APC).^** 



of an ideal Cu02 plane. 

The Knight shift K generally comprises the tempera- 
ture (r)-dependent spin part Ks and the T-independent 
orbital part K^rh as follows: 



(1) 



where a is the direction of B. The spin part of the Knight 
shift for the B \\ a6-plane {K^^) is obtained by subtract- 
ing i^OTbi which is approximately 0.23(±0.02)% assum- 
ing K^^ « in the T=0 limit. For multilayered cuprates 
that exhibit an AFM order, Kf{T) at nonmagnetic OP 
shows an upturn below Tn due to the transferred hy- 
perfine magnetic field arising from AFM IP, as shown in 
the lower panel of Fig. 6: K^"^^ ~0.23(±0.02)% was as- 
sumed to be the same as that of multilayered cuprates in 
which all Cu02 planes are in a paramagnetic state. Note 
that -ft'orV^ ^''^ "^o* different among high-Tc cuprates 
whether at IP or □P.^^-^'^'^g-Ti 



2.3 Hyperfine magnetic field in Cu02 plane 

According to the Mila-Rice Hamiltonian,^"^ the spin 
Knight shift of Cu in the Cu02 plane is expressed as 



(2) 



where Aa and B are the on-site and supertransferred 
hyperfine fields of Cu, respectively. Here, Aa consists 
of contributions induced by on-site Cu 3d.j;2_y2 spins — 
anisotropic dipole, spin-orbit, and isotropic core polar- 
ization, and the B term originates from the isotropic 4s 
spin polarization produced by four neighboring Cu spins 
through the Cu{3dj.2 _y2)-0(2pa)-Cu(4:s) hybridization. 
Since the spin susceptibility Xs{T) is assumed to be 
isotropic, the anisotropy A of Kf{T) is given by 

^ _ KliT) A, + AB 



Kf{T) Aab + 'iB' 



(3) 



The on-site hyperfine fields Aab ~ 3.7 T//iB and Ac ~ 
— 17 T//iB^^ ''^ are assumed as material- independent in 
hole-doped high-Tc cuprates. In multilayered compounds 
[Hgl245(OPT)ttl], B(IP) « 6.1 T/^b and B{OV) « 7.4 
T//XB are estimated, which are larger than B ~ 4 T/^b 
for LSC0,f^8 YBCOt,'^!-'^^''^* and YBa2Cu408(Y1248)^9 
compounds. In the AFM ordered state where spins 
align antiferromagnetically among nearest-neighbor Cu 
sites in a two-dimensional lattice of the Cu02 plane, 
the internal field at the Cu site is generally given by 
Bi^t^\Aab - 4B|A/afm from the AFM moment (Mafm) 
through those hyperfine interactions. 

2-4 Evaluation of planar Cu02 hole density 

In hole-doped high-Tc cuprates, the hole density p 
in Cu02 planes has been determined by various meth- 
ods: indirect chemical methods like solid solutions, ^""^"^ 
bond valence sums determined from structural bond 
lengths,®'*"*^ or the Fermi surface topology.*^ Besides, 
the thermoelectric power is a universal function of p,^^' 
and the phase diagram for hole-doped cuprates is well 
described by Tc(p)=Tc'"^'^[l - 82.6(p - 0.16)2],90 which 
are applicable to the estimation of p when no suitable 
structural data are available. These methods are, how- 
ever, inapplicable to multilayered cuprates composed of 
more than two inequivalent Cu02 planes in a unit cell, 
because these methods evaluate the total hole density not 
each hole density inherent in Cu02 planes. 

In NMR experiments, the spin part of the Knight shift 
at room temperature {K^^iWT)) increases with p from 
the underdoped region to the overdoped region in hole- 
doped cuprates,^^'^^'^^"^^ suggesting that Kf{KY) is 
available to determine planar Cu02 hole densities (p). 
The linear equation p'=0.502-ft:,''''(RT)-f0.0462 has been 
reported, as indicated by the dotted line (Fl) in Fig. 
5,35.37 .^y}^g].g p' ig derived from the NQR frequencies 
of Cu and O in Cu02 planes. However, when noting 
that i^s(RT) for optimally doped cuprates is empirically 
0.35 - 0.39%, an optimal doping level was evaluated to 
be p'=0.22-0.24,39.43,45 ^^ich is relatively larger than 
the widely accepted optimal doping level in hole-doped 
cuprates, i.e., p ~ 0.16, as shown in the figure. This incon- 
sistency is, in part, due to the calculation that connects 
NQR frequency to hole density. 
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Fig. 5. (Color online) Plots of K^''{KT)s for 0212F,'"' 
Bi2Sr2CaCu208 (Bi2212),'J3 Tl2Ba2Cu06+5 (T12201),S* and 
TlSr2CaCu207_5 (T11212),95 as functions of p evaluated us- 
ing Tc=T™*"'[l - 82.6(p - 0.16)2].''° Here, the above cuprates 
are selected because their local structures are homologous to 
those of the multilayered 02(n-l)nF and M12{n-l)n scries, 
which guarantees that the plot of K°^(RT) vs p is material- 
independent among the homologous compounds. The relation- 
ship F2: (p=0.492Kg''''(RT)-0.023) (solid line) between K^''{RT) 
and p obtained by fitting with the data for K^''{KT)< 0.5% 
allows us to separately estimate p inherent in IP and OP in 
multilayered cuprates. More details of the validity of this re- 
lation are described in the literature.^® Note that the relation 
Fl:p'=0.502ii:s''''(RT)-|-0.0462 used in the previous studies had 
overestimated the hole density by 0.06~0.07.'^^"^^ [cited from 
ref.49] 

Recently, we have investigated the bilayered 
(n=2) apical-fluorine compounds Ba2CaCu204(F,0)2 
(n=2:0212F) over a wide carrier density range, ^® 
which provides an opportunity to reexamine p by 
comparing it with those well established in other 
bilayered compounds, exhibiting the maximum Tc 
at p(r™^'')~ 0.16. Figure 5 shows plots of Kf(KT) 
for 7i==2:0212F, BisSrsCaCuzOg (n=2:Bi2212),93 
TlsBaaCuOe+i (n=l:T12201),9** and TlSrsCaCuzOT^^ 
(7i=2:T11212)^^ as a function of p evaluated using 
r^(p)=r^niax[;L _ 82.6(p - 0.16)2],9" indicating that 
Kf{RT) monotonically increases with p from the 
underdoped region to the overdoped region. Here, the 
above cuprates are selected because their local structures 
are homologous to those of the multilayered 02(n-l)nF 
and M12{n-l)n series, which guarantees that the plot 
of K^^{RT) vs p is material- independent among these 
homologous compounds. More details of the validity of 
this relation have been described in the literature.'*^ This 
renewed relation, based on the Knight shift, enables 
us to separately estimate p for each Cu02 plane in 
multilayered cuprates. In this review, we use the linear 
function of p=0.492ii:,"''(RT)-0.023, as shown by the 
solid line (F2) in Fig. 5, which was obtained by fitting 
with the data for Kf{RT)< 0.5%, since the Kf(RT) 
of the samples used in this study is less than 0.5%. Note 
that the relation (Fl) used in the previous studies had 



overestimated the hole density by 0.06^0.07.'^^""' 

According to eq. (2), the p dependence of Kf^(RT) 
is derived from those of the B term in the hyperfine 
coupling constant and Xs(RT). In multilayered cuprates, 
the B term increases moderately with p, showing a 
steep increase at p=0.18^0.20.^^ The B term arises 
from Cu{3dx2^y2)-0{2pa)-Cu{4:s) covalent bonds with 
four nearest-neighbor Cu sites; therefore, the large B sug- 
gests the strong hybridization between the Cu{3dx2-y2) 
and 0{2pa) orbits. This result is consistent with the fact 
that a metallic state is more stabilized in an overdoped 
regime. For layers with tetragonal symmetry in cuprates 
homologous to the multilayered series, the p-dependent 
B terms are material-independent j"^^ being larger than 
the B terms for LSCO, YBCOe+x, and Y1248, - 4 
T//iB.^*' This difference is also seen in the variation in 
the nuclear quadrupole frequency ^^vq. ^^vq increases 
with p for all the materials,^^'^^'^^ while, for a given p, 
the absolute values of what LSCO and YBCOg-i-a; are 
about 2 to 3 times larger than those for others. ^^'^Q 
depends on the hole number in a Cu{3dx2-y2) orbit 
and on rip in a 0{2p(7) orbit. Therefore, it is expected 
that Ud and Up in LSCO and YBCOg+x will be different 
from those in the other compounds treated here, even if 
they have the same p {p^nci+2np—l) . Actually, it has 
been reported that Ud is large in LSCO and YBCOe+x, 
which is the reason for the large i/q.^^'^'^ In this context, 
the reason why the B terms in LSCO, YBCOe+x, and 
Y1248 deviate from those in other materials is the dif- 
ferent partitions of holes in the Cn{3dx2-y2) and 0(2pcr) 
orbits, which is probably related to the crystal structures; 
LSCO, YBCOg-i-a;, and Y1248 have orthorhombic crystal 
structures in the superconducting region, while 0212F, 
Bi2212, T12201, and T11212 have tetragonal ones. Thus, 
we conclude that the p dependence of B holds in Cu02 
planes with tetragonal symmetry homologous to the mul- 
tilayered series, which guarantees the estimation of p for 
OP and IP independently in multilayered cuprates using 
the renewed relation between K^''{KT) and p. 

3. Results 

3.1 Five-layered (n—5) compounds 

Figures 6(a) - 6(g) respectively show the layer- 
dependent physical properties unraveled by site-selective 
NMR studies of the following n—5 compounds: 
Hgl245(UD)ttl (Tc=72 K),39'5i Hgl245(UD)p (Tc=82 
j^) 46,51 Hgl245(UD)p (Tc=92 K),*^^^! Hgl245(OPT)ttl 
(Tc=108 K),38.43 Hgl245(OPT)tt2 (Tc^llO K),^^ 
T11245(OVD) (Tc^lOO K),38.43 and Cul245(OVD) 
(Tc=90 K).^'' The temperature (T) dependences of 
K'f'{T) are shown in the lower panel of the figures. The 
''(RT)s at room temperature decrease with decreasing 
p, which enables us to evaluate p at each Cu02 plane 
using the relationship of ^"''(RT) vs p.^'^ 

3.1.1 Superconducting characteristics 

In multilayered cuprates, TcS inherent in OP and IP 
are estimated from the T dependences of K^'^ and its 
T derivatives,'^^"'^'''^'^ which are shown in Figs. 7(a) and 
7(b). Figure 7(a) shows that a distinct peak in the T 
derivatives of K^^ at OP coincides with the bulk Tc= 110 
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Fig. 6. (Color online) Illustration of layer-dependent physical properties for n=5 compounds: (a) Hgl245(UD)Hl,^^' (b) 
Hgl245(UD)tt2,46. 51(c) Hgl245(UD)t)3,*®'5i (d) Hgl245(OPT)t)l,3*'*3 (e) Hgl245(OPT)|t2,43 (f) T11245(OVD),=*8'43 and (g) 
Cul245(OVD).'^^ The middle panels present tables of the hole densities of p(IP) and p(OP), and the AFM ordered moments of 
A/afm(IP) and AfAFM(OP). The lower panels show the T dependences of K^^{T)s, which enable us to separately estimate ps for IP 
and OP, and to probe the onset of AFM and HTSC at IP and OP. 



K in optimally doped Hgl245(OPT)jJ2, demonstrating 
that the SC is driven by OPs with an optimum p, while 
another peak at r=85 K in the T derivatives of K^^ at 
IP is assigned to inherent to IP with a hole density 
smaller than that at OP.''"^ 

Figure 7(b) shows that a distinct peak in the T deriva- 
tives of Kf at IP coincides with the bulk Tc = 90 K in 
overdoped Cul245(OVD), demonstrating that the SC is 
driven by IPs, while a distinct peak at T = 65 K in the 
T derivatives of K^'' at OP was assigned to inher- 
ent in overdoped OP."^^"'^^ When noting that OP and IP 
are alternatively stacked along the c-axis in multilayered 
cuprates, the SC gaps at IP (OP) for Hgl245(OPT)tJ2 
[Cul245(OVD)] with lower than bulk are antici- 
pated to develop linearly between Tc and T^, as shown in 
Fig. 7(c) [Fig. 7(d)], respectively, owing to the proximity 
effect. Likewise, it was widely established that Tc in mul- 
tilayered cuprates is determined by the carrier density at 
each Cu02 plane,'^^"'^^''*'^''*^'^^ and it was proposed that 
such two-gap SC causes anomalous behaviors. 

3.1.2 Estimation of AFM moments ~ zero-field Cu- 
NMR/NQR studies 
The observation of a zero-field (ZF) Cu-NMR spec- 
trum at low T allows us to estimate the AFM ordered 
moment Mafm- In general, the Hamiltonian for Cu nu- 
clear spins (/ = 3/2) in crystal lattices with an axial 
symmetry is described by the Zeeman interaction due to 
the magnetic field B, Hz, and the nuclear-quadrupole 



interaction T-Lq as 
H = -Hz + Hq 

= -^^;,J.B + _|^(34_/(/+i))(4) 

where is the Cu nuclear gyromagnetic ratio, eQ is 
the nuclear quadrupole moment, and eq is the electric 
field gradient at the Cu nuclear site. In T-Lq, the nu- 
clear quadrupole resonance (NQR) frequency is defined 
as lyq = e^qQ/2h. In nonmagnetic substances, the NQR 
spectrum originates from the second term in eq. (4) at 
a zero external field (-Bcxt=0 T). On the other hand, in 
magnetically ordered substances, the ZF-NMR spectrum 
is observed owing to the internal magnetic field i?int at 
Cu sites in eq. (4) despite i?oxt=0 T. Figure 8(a) shows 
a typical Cu-NQR spectrum of n=4:Hgl234(OPT) with 
Tc=123 K at 1.5 K.'^'^ The respective ^^Cu-NQR frequen- 
cies i^^vg) at IP and OP are 9.6 and 17.8 MHz, which 
are comparable to the typical ^^i^q values at IP and OP 
in a paramagnetic regime. ^^'^^ 

Figures 8(b)-8(g) respectively show the Cu-NQR/ZF- 
NMR spectra at 1.5 K for T11245(OVD)38.43 
Hgl245(OPT)p,43 Hgl245(OPT)tJl,38.43 
Hgl245(UD)p,5i Hgl245(UD)tt2,5i and 

IIgl245(UD)ttl.^^'^i In these spectra, no NQR spectrum 
at IP is observed, pointing to an onset of AFM order 
at IP, whereas the NQR spectrum at OP is observed 
for (b) T11245(OVD),38.43 (g) Hgl245(OPT)p,43 
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Fig. 7. (Color online) T dependences of K^^s (solid circles) 
and its T derivatives (empty circles) at OP and IP of (a) 
Hgl245(OPT)t)2 with rc=110 K and (b) Cul245(OVD) with 
Tc=90 K. Owing to a large imbalance in ps , TcS inherent in 
OP and IP are assigned from a peak in the T derivatives of K^^ 
(see text), (a) A distinct peak in the T derivatives of K^^ at OP 
coincides with the bulk Tc= 110 K, demonstrating that the SC 
is driven by OPs. Another peak at T=85 K in the T derivatives 
of K^'' at IP was assigned to inherent in IP.''^ (b) A distinct 
peak in the T derivatives of K^^ at IP coincides with the bulk 
Tc= 90 K, demonstrating that the SC is driven by IPs. A distinct 
peak at T=65 K in the T derivatives of K^'' at OP was assigned 
to inherent in overdoped OP.^^"^'^ (c,d) The SC gaps at IP 
(OP) for Hgl245(OPT)tt2 (Cul245(OVD)) with T^' lower than 
bulk Tc are anticipated to develop linearly between Tc and 
owing to the proximity effect.^® [cited from refs.^^'*'^] 



(d) Hgl245(OPT)ttl,38'43 (e) Hgl245(UD)tt3, and (f) 
Hgl245(UD)tJ2. 

These spectra at IP, that are observed in the range 
of 20^^50 MHz, are reproduced by assuming the in- 
ternal field Sint at IP, which is generally given by 
Ant=|^hf|MAFM=|^ah - 4S|A/afm. Here, Aab « 3.7 
T/^B, B{IF) w 6.1 and B(OP) « 7.4 T/^b are 

assumed in multilayered compounds,'^® Mafm is a spon- 
taneous AFM moment at Cu sites. Using these values, 
Mafm(IP)s at IP are estimated to be in the range of 
0.1^0.18 fiB at r=1.5 K, which are listed up in the 
middle panel of Fig. 6. Note that the Cu-NMR spectra 
in undoped Mott insulators are observed in the higher- 
frequency ranges, such as 75'^110 MHz for La2Cu04^''^ 
and YBC06^°^'^°^ or 125-150 MHz for CaCu02(IL) 
(see Fig.l9(f)) from the MAFM=0.5-0.7/iB.^°^' Thus, 
the mobile holes existing at IP uniformly reduce to 
Mafm(IP)=0.1~0.18 hb, indicating that a static AFM 
metallic state is realized at IP uniformly. Here, note that 
the presence of the spin-glass phase is excluded because 
the possible distribution in A/afm(IP) is less than ± 
0.02/ZB. 

In the most underdoped Hgl245(UD)jJl, the spec- 
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Fig. 8. (Color online) Cu-NQR/zero-field NMR spectra at 
1.5 K for (b) T11245(OVD),38>« (c) Hgl245(OPT)tt2,« 
(d) Hgl245(OPT)ttl,38,43 (e) Hgl245(UD)tt3,5i (f) 
Hgl245(UD)tt2,5l and (g) Hgl245(UD)ttl,33' along 
with (a) n=4:Hgl234(OPT) with Tc=123 K in the param- 
agnetic state. Dotted lines represent the NQR frequencies 
63vq(IP)=8.4 and ^'■^vq{OV)=\Io MHz for Hgl245(OPT)ttl.='8 
The solid bars represent resonance frequencies and intensities 
for two components of isotopes, ^^Cu and ^^Cu, which were 
calculated based on eq. (4). [cited from refs.^*' ^'^J 



trum at OP is observed at approximately 30 MHz (see 
Fig. 8(g)), indicating that Mafm(OP)'- 0.092/iB even at 
OP that is responsible for the onset of HTSC with T^— 
72 K. Both the phase separation and the spin-glass phase 
are excluded since no NQR spectra are observed. These 
facts indicate that both AFM and HTSC uniformly coex- 
ist at the microscopic level. '^^ Here, note that the broad- 
ening of these spectra at OP and IP points to a possi- 
ble distribution of Mafm (OP) and Mafm (IP) with ap- 
proximately ± 0.02 \iB over the samples. This may be 
partly because some subtle inhomogeneity in carrier den- 
sity takes place owing to deoxidization process. Note that 
Mafm and Tc at OP in Hgl245(UD)ttl are comparable 
to those at IP in Hgl245(OPT)t|2, since the p values are 
almost the same for these layers. This means that Mafm 
and Tc are primarily determined by p at each layer not by 
the electronic states of adjacent layers, whereas Tn de- 
pends on the number n of Cu02 layers or on interlayer 
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magnetic coupling, as discussed in §4.1 and §4.2. 

3.1.3 Determination of Neel temperature Tn for AFM 
order 

The Neel temperature Tn is determined by the mea- 
surement of nuclear-spin-relaxation rate 1/Ti, which ex- 
hibits a peak at Tn- Generally, l/Ti is described as 

1 _ ^ , ^ |2 lm[x(g:^o)] 

where Aq is a wave-vector (q)-dependent hyperfine- 
coupling constant, x(<7,a;) is the dynamical spin suscep- 
tibility, and ujq is the NMR frequency. The T depen- 
dence of 1/Ti exhibits a peak at Tn because low-energy 
spectral weight in x{q = Qi'^) is strongly enhanced at 
Wo ~ in association with the divergence in magnetic 
correlation length at T Tn- Here, Q is an AFM wave 
vector (tt, tt). In the case of Hgl245(OPT)ttl, for in- 
stance, the AFM order at IP was detected by ^^Cu-NMR 
l/Ti at OP that shows a peak at rN=60 K, as shown in 
Fig. 9(b). 3* Likewise, the T^s at IP were 55 and 45 K 
for Hgl245(OPT)tt243 and T11245(OVD),38'42,43 respec- 
tively. Note that the AFM order below Tn = 60 K in 
Hgl245(OPT)t|l was also confirmed by the /iSR measure- 
ment. Note that the onset of AFM order at IP was 
corroborated by an upturn in K'f'{T) at OP upon cool- 
ing below the Tns in Hgl245(OPT)jJl, Hgl245(OPT)jl2, 
and T11245(OVD), as marked by upward-arrows in the 
lower panel of Fig. 6. As a result, the Tns for the 
more underdoped compounds such as Hgl245(UD)t)l, 
Hgl245(UD)tJ2, and Hgl245(UD)p were evaluated from 
the temperature below which the Kf'{T) at OP reveals 
an upturn. As expected, as p decreases, the Tn at IP in- 
creases from Tn -110, 150, and 180 K for Hgl245(UD)p, 
Hgl245(UD)tt2,46^5i and Hgl245(UD)ttl, respectively.^^ 

The T variation of Mafm below Tn was indirectly 
probed on the basis of the T dependence of the inter- 
nal field at nuclear sites in charge reservoir layers. Fig- 
ure 9(a) shows the T dependence of the ^^^Hg-NMR 
spectrum for the B _L c-axis and its full-width at half- 
maximum (FWHM) of Hgl245(OPT)ttl with Tc=108 K 
and Tn=60 K. The FWHM increases rapidly below Tn= 
60 K, probing the development of the internal field at the 
iQ^Hg site induced by the onset of Mafm (IP)- The T de- 
pendence of Mafm (IP) was close to the theoretical pre- 
diction for weak itinerant AFM metals, as indicated 
by the solid line in Fig. 9(a). Note that the FWHM of 
the ^^^Hg-NMR spectrum is increased markedly below 
Ta= 25 K, which has not been identified yet.'^* 

Recently, in the apical-F multilayered compound 
n=5:0245F with Tc=52 K, which is more underdoped 
than n=5:Hgl245(UD)Ul with Tc=72 K, clear evidence 
of the T evolution of Mafm bas been presented be- 
low Tn=175 K, along with an SC diamagnetic shift 
below Tc=52 K through ^^F-NMR studies.-^^ The ZF 
Cu-NMR study revealed that Mafm(IP)=0.20 /ib and 
Mafm(OP)=0.14 /iB at 1.5 K. As shown in Fig. 10, 
the internal field |i3f„t(F)| at the apical-F site for the 
B II c-axis, which was evaluated from the splitting of the 
^^F-NMR spectra, increases significantly below Tn = 175 
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r(K) 

Fig. 9. (Color online) (a) T dependences of ^^^Hg-NMR spec- 
trum and its full-width at half-maximum (FWHM) when the 
B ± c-axis for n=5:Hgl245(OPT)t)l. The solid line indicates 
the fitting of FWHM to the relation MA¥yi{T) oc (1 - 
(T/Tjv)^'''^)^/'^ , which is the theoretical prediction for weak itin- 
erant AFM metals. (b) T dependences of short (triangles) and 
long (squares) components in ®^Cu-NMR 1/Ti at OP. A peak 
in 1/Ti points to Tn=60 K, below which the FWHM at the Hg 
site increases owing to the development of A/afmCIP)-^* [cited 
from ref.'^*] 

K,"° at which the ^^F-NMR 1/Ti also exhibits a peak. 
The T dependence of |i3f„j(F)| at the apical-F site for 
B II c-axis was roughly reproduced down to Tc— 52 K by 
either Sf„t(F)cx Mafm(T)=Mafm(0)(1 - (T/Tw)3/2)i/2 
(solid line)i°9 or oc (I - {T/Tn)Y/^ (dotted line) used 
for slightly doped LSCO compounds. These results 
convinced us of the three-dimensional long-range AFM 
order below Tn=175 K in n=5:0245F with the SC tran- 
sition below Tc=52 K. Here, we note that B^^^{¥) shows 
an additional increase as T decreases below Tc; the in- 
crease in B^^^^iY) below Tc implies that the onset of an SC 
order parameter is actually coupled with Mafm (OP) in 
the HTSC-AFM coexisting state."*^ This result gives con- 
vincing evidence that the HTSC below Tc=52 K emerges 
on the background of the AFM order taking place below 
Tn=175 K; hence, both coexist. 

3.1.4 Pseudogap behavior in n—5 compounds 

Since the discovery of high-Tc cuprates, anomalous 
normal states exhibiting a pseudogap behavior have been 
one of the most important subjects in the research on 
HTSC. Initially, a gaplike behavior was reported as a 
gradual suppression of 1/TiT below T*,^^^ which has 
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Fig. 10. (Color online) (a) T dependence of the internal field 
at apical F site (|Bfnt(F)|) for n=5:0245F with Tc=52 K and 
Tn=175 |Bf^t(F)| was estimated from the splitting of the 

l^F-NMR spectra in the B \\ c-axis. (b) Tn=175 K was de- 
termined by a peak of ^^F-NMR l/Ti. The sohd and dot- 
ted lines represent B9^^{F)oc (1 - (T/Tiv)^/^)^!/^' and oc 
(1 -T/Tiv)(i/2) ^111 respectively, [cited from ref.^S] 
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been called the spin gap. Neutron scattering experi- 
ments in underdoped regions also showed that spin ex- 
citations at low energies are suppressed in the nor- 
mal state. ^^■^ Then, angle-resolved-photoemission spec- 
troscopy (ARPES) experiments directly identified the ex- 
istence of an energy gap in electronic spectra even above 
114,115 r^Yns gap observed in ARPES below T* turned 
out to have the same angular dependence as the d-wave 
SC gap in the Brillouin zone.^^^'^^^ After these exper- 
imental observations of single-particle spectra, the gap 
has been called the pseudogap. 

Figure 11 indicates the T dependence of ^'^ Cu- 
NMR 1/TiT at IP and OP of the n=5 compounds 
for the B \\ c-axis.37.38,42,43 r^Yie l/TiT at OP for 
Hgl245(OPT)ttl starts to decrease upon cooling below 
r* = 160 As shown in Fig. 12, T* decreases as p in- 
creases as is also observed in single- and bilayered com- 
pounds.^'^' ^^^' However, note that when the under- 
doped IPs of Hgl245(OPT)ttl and T11245(OVD) show 
an AFM order below Tn=60 K and 45 K, respectively, 
l/TiT exhibits no indication of a pseudogap behavior 
above 140 K. These results reveal that low-energy spec- 
tral weights in at IP of Hgl245(OPT)t|l and 
T11245(OVD) are critically enhanced around a; —> to- 
ward the AFM order. As a result, the NMR spectra at IP 
are lost below 140 K owing to the extremely short nuclear 
relaxation times. By contrast, l/TiT at the slightly over- 
doped IP for Cul245(OVD) shows the pseudogap state 
below T*=145 K.^'^'^^'''^ Eventually, we highlight that a 
spin-gap collapses in the underdoped region where the 
AFM order emerges. Recently, a similar behavior on the 
p dependence of T* has been observed in the n=3 com- 



Fig. 11. (Color online) T dependence of ^^SCu-NMR 1/TiT at IP 
and OP in n=5 compounds.''^' The data are presented 
from the most underdoped IP (top) to the heavily overdoped OP 
(bottom). The T* of n=5 compounds increases as p decreases, 
but the underdoped IP, which shows the AFM order below 60 K 
(Hgl245(OPT)ttl) and 45 K (T11245(OVD)), shows no indication 
of a pseudogap above 140 K. This suggests that the spin-gap 
collapses in the AFM-HTSC mixed state for p < pc('«). [cited 
from refs.37.38.42,43] 

pounds.^" This pseudogap behavior is also extensively 
discussed in § 4.4. 

3.1.5 Phase diagram in n=5 compounds 

Figure 12 shows a novel phase diagram for n=b com- 
pounds. Here, Tn and Tc are plotted as functions of p, 
which is estimated from the relationship of K^^{KT) vs 
p discussed in § 2.4. The characteristic features are sum- 
marized as follows: (i) The AFM metallic phase (AFMM) 
is robust up to p ^0.1 and uniformly coexists with the 
HTSC state up to a quantum critical density pc{n=b) 
0.1 at which the AFM order collapses, (ii) Tc has a peak 
at p ~0.16 after the AFM order collapses. These find- 
ings suggest the intimate relationship between HTSC 
and AFM. (in) The p dependence of T* below which 
l/TiT starts to decrease points to that a spin-gap col- 
lapses in the AFM-HTSC mixed state. These features 
of a phase diagram differ significantly from the well- 
established phase diagram of LSCO,'^'' in which the AFM 
and HTSC phases are separated by the spin-glass phase 
in association with the Anderson localization mechanism 
(see Fig. 20). 
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Fig. 12. (Color online) Phase diagram of AFM and HTSC at ho- 
mogeneously doped Cu02 plane in n=5 compounds. Tn, Tc, and 
T* are plotted against the hole density p.^^ The solid and empty 
circles correspond to the data for IP and OP, respectively. The 
AFM metallic phase is robust and uniformly coexists with the 
HTSC state up to a quantum critical density Pc(5) ~ 0.1 at which 
AFM order collapses. Tc exhibits a maximum at p(T^'^^) ~ 0.16. 
Here, PM and PC denote the paramagnetic phase and pseudogap 
state, respectively. The star denotes the pseudogap temperature 
T*. [cited from refs.39.«.49] 



3.2 Phase diagram in n=4 compounds 

Since it was difficult to change carrier density widely 
in n=4:Hgl234,47 we used Ba2Ca3Cu408(FyOi_j,)2 
(n=4:0234F)^^ as n—A compounds. The substitution of 
for at the apical site increases p and Tc from 55 
K up to 102 K.^'S'S^ As shown in Figs. 13(a)-13(d), sys- 
tematic ^^Cu- and ^^F-NMR studies have revealed SC 
and AFM characteristics for n=4:0234F, which are de- 
noted as 0234F(ttl), 0234F(tt2), 0234F(|13), and 0234F(tt4) 
with nominal contents y= 0.6, 0.7, 0.8 and 1.0, respec- 
tively. The measurements of K^^ shown in Figs. 13(e-h) 
enable us to estimate p at IP and OP for the samples. 
In fact, p increases progressively with decreasing y. 

Figures 14(a)-14(d) show the Cu-NQR/ZF-NMR spec- 
tra at r=1.5 K for n=4:0234F: (a) tJl(y --0.6), (b) 
p(y -0.7), (c) tt3(y -0.8), and (d) ^(y -1.0). For ^1, 
two NQR spectra revealed that ^^iyQ{lP)=9.7 MHz and 
^'^z/q(0P)=15 MHz, which are comparable to those for 
other paramagnetic multilayered cuprates.**^'^^'^^ For p, 
the NQR spectrum at OP is observed at ^^i/q(0P)=15 
MHz. By contrast, it is not expected that the NQR spec- 
tral intensity at IP(i) observed at 9.1 MHz will be signif- 
icantly smaller than the intensity of the ZF-NMR spec- 
trum at IP(ii). The latter spectrum probes an internal 
field i3int=1.5 T in association with the onset of AFM 
order with Mafm(IP)~ 0.08 fiB- These results point to 
that the IPs for p undergo a phase separation into the 
paramagnetic and AFM phases owing to closeness to Pc, 
at which an AFM order collapses. 

For (13, the observation of the NQR spectrum at OP 
indicates that a spontaneous AFM moment is absent, 
whereas the NMR spectrum at approximately 28 MHz 
at IP probes Bint ^ 2.4 T and hence Mafm(IP)-- 0.12 
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Fig. 13. (Color online) Illustration of layer-dependent physical 
properties for n=4:0234F; (a) J4(i/ ~1.0), (b) tt3(?/ ~0.8), (c) 
i2{y ~0.7), and (d) il(y ~0.6). The middle panels present tables 
of the hole densities of p(IP) and p(OP), and the AFM ordered 
moments of Mafm(IP) and Mafm(OP) (see text). The lower 
panels show the T dependences of K^^{T)s, which enable us to 
separately estimate ps for IP and OP, and to probe the onset of 
AFM and HTSC at IP and OP. [cited from refs.''5.49] 



^B- Figure 14(d) shows that the Cu-ZF-NMR spectra at 
IP and OP are observed at approximately 45 and 30 MHz 
for t|4 with Tc=55 K, which allows us to estimate -Bint ^ 
3.8 T and 2.7 T and hence Mafm(IP)'- 0.18/ib and 
-/^afm(OP)— 0.11/xb, respectively. Note that no trace 
of NQR spectra excludes the possibility of phase separa- 
tion into the paramagnetic and AFM phases. Therefore, 
the OP, which is mainly responsible for the HTSC with 
Tc= 55 K, manifests the AFM order, indicating that the 
uniform mixing of AFM with the spontaneous moment 
Mafm=0.11/ib and HTSC at Tc= 55 K occurs in the OP. 
Taking into account of the ARPES experiment on this 
compound (P) that observed SC gaps on Fermi sheets 
of both IP and OP,"^ we deduce that the AFM order 
coexists uniformly with HTSC at both IP and OP for 

tt4.53 

The TnS of tt2, p, and tJ4 are determined by the ^^F- 
NMR l/Ti measurements at the B \\ c-axis, as presented 
for all the samples in Figs. 14(e-h). In the present case, 
^^(l/Ti) is dominated by magnetic fluctuations.^^ Fig- 
ures 14(f) and 14(g) show the T dependences of the 
^^(l/ri)s of p and p, which exhibit peaks at approx- 
imately 30 and 50 K, respectively. These peaks are as- 
sociated with the onset of the AFM order at Tn=30 K 
at IP(ii) for p and Tn=50 K at IP for p, accompa- 
nying MAFM(IP(ii))--0.08 fiB and Mafm(IP)--0.12 ^b, 
respectively. It is unexpected for the 1/Ti of P to show 
two peaks at Tn — 80 K and T-^ — 30 K, as shown 
in Fig. 14(h). Since p{lP) < p{OF) and Mafm(IP) > 
A^afm(OP), it is likely that TnS inherent to IP become 
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Fig. 14. (Color online) Cu-NQR/ZFNMR spectra at Boxt=0 T 
and T=1.5 K for n=4:0234F: (a) ftl, (b) p, (c) P, and (d) p. 
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and (h) p at /=174.2 MHz and B || c-axis. The AFM moments 
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text), [cited from ref.^^] 
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Fig. 15. (Color online) Phase diagram of AFM and HTSC at 
homogeneously doped Cu02 plane in n=4 compounds. Tn and 
Tc are plotted against the hole density p. The solid and empty 
circles correspond to the data for IP and OP, respectively. Note 
that the uniform mixing of AFM with Tn=30 K and HTSC with 
Tc=55 K takes place at OP of p. A critical hole density for the 
AFM order in n=4 compounds is extrapolated to Pc(4) ~ 0.08, 
which is lower than that for the n=5 compounds. The pseudogap 
phase of the n=4 compounds has not been determined yet. [cited 
from refs.*^'*^l 



deviation from the nominal apical-fluorine F ^ content 
results in the doping of hole carriers into both OP and 
IP. 



larger than that at OP. The AFM order inherent to OP 
may develop below ~ 30 K. Here, note that T^^ ^ 
30 K at OP for |14 is comparable to Tn 30 K at IP 
(ii) for p because their ps are almost the same for both 
compounds. 

Figure 15 shows the phase diagram of AFM and HTSC 
for n=4:0234F. Here, and Tn at OP and IP are plot- 
ted against p.^^ This phase diagram resembles that for 
the n=5 compounds. In particular, the uniformly mixed 
phase of AFM with Tn=30 K and HTSC with Tc=55 
K was observed at the OP for jj4, demonstrating that it 
is the universal phenomenon inherent in a single Cu02 
plane in the underdoped regime. In the phase diagram of 
the n=4 compounds presented in Fig. 15, Pc(4) at which 
the AFM order collapses is extrapolated to 0.08, which 
is smaller than Pc(5) ~0.1 for the n=5 compounds. As 
discussed in § 4.1, this is because the interlayer magnetic 
coupling becomes weaker as the number n of Cu02 layers 
decreases. 

Here, we comment on 0234F(tj4) that attracted much 
attention at the initial stage of studies, since the 
self-doping mechanism had been proposed: charge carri- 
ers are transferred between IP and OP in a unit cell. 
Here a formal Cu valence was assumed to be just 2+ in 
an ideal case of a nominal fluorine content (F2), where 
the apical sites are fully occupied by F~^. However, ex- 
tensive investigations on 02(n-l)nF with n=2, 3, and 4 
have ruled out the possibility of a self-doping mechanism 
in these compounds.'*'* The point was that an inevitable 



3.3 Phase diagram ofn—3 compounds 

The 63Cu- and ^^F-NMR studies of three-layered 
Ba2Ca2Cu306(Fj^Oi_j,)2 compounds with a nominal F 
content y denoted as rt=3:0223F have revealed the SC 
and AFM characteristics shown in Figs. 16(a)-16(d).'**'^° 
The ps at IP and OP determined by the relationship 
of KfiRT) vs p are listed in the middle panel of 
Fig. 16. Figure 16(e) shows the T dependences of K^'^ 
and dKf/dT for p. Besides the bulk Tc=76 K, in- 
herent in IP is tentatively deduced as Tp'=60 K from a 
secondary peak in the T dependence in dK^^ /dT. Fig- 
ure 16(f) shows the Cu-NQR/ZFNMR spectrum of U4, 
along with that of jjl; the spectra are observed at each 
NQR frequency, revealing no AFM order at either layer 
of jjl. In tt4, the spectrum observed at 13.6 MHz cor- 
responds to the vq at OP, which is slightly lower than 
that at OP of jjl because of the lower doping level. 
On the other hand, the spectrum of IP is observed at a 
frequency much higher than the NQR frequency, which 
probes Bint ^ 2.4 T and hence Mafm(IP)-^ 0.12 ^b- We 
note that a phase separation into the AFM and param- 
agnetic phases is excluded in the IP of P, because no 
paramagnetic NQR spectrum for IP was observed. The 
Tn=23 K at IP was determined from the peak in *^F- 
NMR l/Ti (see ref.*"). Details of systematic NMR exper- 
iments on jjl, p, and p were published elsewhere. As a 
consequence, we have unraveled a phase diagram of the 
n=3:0223F compounds,^° as shown in Fig. 17. The AFM 
order emerges at IPs of p and P with p < 0.075, suggest- 
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Fig. 17. (Color online) Phase diagram of n=3 compounds. The 
solid and empty circles correspond to the data for IP and OP, 
respectively. The critical hole density of AFM order for the n=3 
compounds is extrapolated to Pc(3) ~ 0.075, which is lower than 
that for the n=4 and n=5 compounds. The filled and open stars 
represent the pseudogap temperature T* determined from the 
peak of 63Cu-(l/TiT) for n=3:0223F50 and M1223,50> 65. "0 ^g. 
spectively. This suggests that the spin gap collapses in the AFM- 
HTSC mixed state for p < pc(3). [cited from refs.'"''50] 



Fig. 16. (Color online) Illustration of layer-dependent physical 
properties for n=3:0223F: (a) tt4(j; ~1.0), (b) P(y ~0.9), (c) 
tt2(y ~0.8), and (d) il{y ~0.6). (e) T dependences of K^'' and 
dKf/dT for p. (f) Cu-NQR/ZFNMR spectrum for (14 at 1.5 K, 
along with that for Jl with no AFM order, revealing Mafm(IP)~ 
0.12 /iB and Mafm(OP)=0 for ((4. The middle panels present 
tables of the hole densities of p(IP) and p(OP), and the AFM 
ordered moments at IPs (see text). [cited from refs.^*' ®'^] 

ing that the critical hole density of AFM order for n—2> 
compounds Pc(3) is close to 0.075, which is lower than 
Pc(4) ~0.08 for the n=4 compounds and Pc(5) ~0.10 for 
the n=5 compounds. These AFM ordered states at IPs 
appear under the background of HTSC with T^— 60 K, 
which indicates that the uniformly mixed state of AFM 
and HTSC emerges universally for a single Cu02 layer 
in the underdoped region, when a magnetic interlayer 
coupling is strong enough to stabilize an AFM order. 

To probe the pseudogap behavior for the n=3 com- 
pounds, the ^'^Cu-NMR 1/TiT measurement has been 
performed in jjl, p., and tt3.^° Fi gure 18 shows the T de- 
pendence of l/TiT for p, which starts to decrease upon 
cooling below T* ^ 150 K for OP. On the other hand, 
the l/TiT at IP in the same sample continues to increase 
upon cooling down to 200 K, and the NMR spectrum at 
IP is lost below 200 K because of extremely short nuclear 
relaxation times due to the critical enhancement of AFM 
spin fluctuations. Taking into account of the doping level 
{p ^ 0.073), unless an AFM order occurs, one could ex- 
pect T* to be around 200 ~ 250 K, which is deduced from 
the case of the single- and bilayered cuprates presented 
in Fig. 24, but no peak of 1/TiT was seen above 200 K. 
Instead, the inset of Fig. 18 indicates that the l/TiT for 
IP at high temperatures follows the Curie- Weiss law as 

l/TiT=C/{T+9) with 9 10 K. This indicates that 

the out-of-plane magnetic interaction, which is responsi- 
ble for the onset of the AFM order at Tn ^ 10 K, causes 



the in-plane AFM correlation to develop further without 
opening a gap at the low-energy spectral weight of AFM 
excitations. Although l/TiT is expected to decrease ow- 
ing to the opening of the SC gap, we consider that the 
low-energy parts in AFM excitations below Tc continue 
to increase at IP down to Tn ~ 10 K, since the peak at 
T ~ 10 K observed in the 1/TiT for OP points to some 
divergence of the 1/TiT at IP towards the AFM order at 
Tn ~ 10 K through the supertransferred hyperfine cou- 
pling between ^^Cu nuclei at OP and the Cu-derived spin 
at IP. When noting that 9 coincides with Tn, as far as 
the low-energy parts in AFM excitations are concerned 
in an energy region smaller than an SC energy gap, we 
suppose that a spin gap may not open between 10 K and 
200 K. 

3.4 AFM and HTSC in n—8 compound 

To gain further insight into the n dependence of Pc{n) 
at which the AFM order collapses, we have investigated 
the magnetic and SC properties of n=8:IIgl278(OPT) 
with rc=103 K}^^ Here, note that the interlayer mag- 
netic coupling is expected to be larger than in the other 
multilayered cuprates. As indicated in Figs. 2(a) and 
2(b), note that the c-axis length increases progressively 
with n, whereas Tc is almost constant for n > 6. The 
sample for NMR study comprises an almost single phase 
of n=8:Hgl278, although it is inevitable that a small 
number of stacking faults along the c-axis are contained 
for layered cuprates with n > 6.^'''^" 

Figures 19(a)-19(c) show the T dependence of 
K^'' at OP of (c) n=8:Hgl278(OPT), together with 
(a) n=4:Hgl234(OPT) and (b) n=5:Hgl245(OPT)tJl, 
which are as-prepared samples in an optimally doped 
regime. X°''(RT) at OP is comparable to those for 
n=4:Hgl234(OPT) and n=5:Hgl245(OPT)ttl, pointing 
to p(OP) ^ 0.14 in these compounds. It is also corrob- 
orated by similar NQR frequencies of OP in these com- 
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Fig. 18. (Color online) T dependences of ''^Cu-NMR 1/TiT at 
OP and IP for {13 at _B ± c and /=174.2 MHz.^'^ 1/TiT at 
IP continues to increase upon cooling towards the AFM order at 
Tn ~ 10 K, whereas 1/T\T at OP exhibits a pseudogap behavior 
with T* '^150 K. The p dependence of T* for ?i=3 compounds is 
shown in Fig. 17. The inset shows a plot of TiT vs T to reveal the 
Curie- Weiss behavior at high temperatures as 1 /TiT=C / {T + d) . 
As for IP, 6=-10 K was obtained, which corresponds to the peak 
of l/TiT at OP, that is, ~10 K at IP.^O [cited from rei.'^°] 



pounds as will be discussed later. Therefore, if all lay- 
ers were in the paramagnetic state, the T dependence 
of Kf should be similar to that of n=4:Hgl234(OPT) 
as indicated by the broken curves in Figs. 19(b) and 
19(c). By contrast, note that the Kf at OP for 
n=5:Hgl245(OPT)ttl deviates from this broken curve 
below Tn=60 K at which the AFM order sets in at 
IPs. Accordingly, the Tn of n=8:Hgl278(OPT) is ten- 
tatively deduced to be about 300 K, since K^'' devi- 
ates from the broken curve below 300 K, as shown in 
Fig. 19(c). Although n=8:Hgl278(OPT) comprises three 
crystallographically inequivalent IPs, no Cu-NMR spec- 
tra at IPs were observed up to 300 K, which suggests 
that AFM order emerges at all IPs. To evaluate AFM 
moments of IPs, we present the Cu-NQR/ZF-NMR spec- 
tra of Hgl278(OPT) at r=1.5 K in Fig. 19(f), which 
are compared with the spectra of (d) n=4:IIgl234(OPT) 
and (e) n=5:Hgl245(OPT)ttl. Since vq{OV)^ 16 MHz at 
OP for n=8:Hgl278(OPT) coincides with those for n=4 
and n=5 compounds, the OP is paramagnetic. The Cu- 
ZF-NMR spectrum at approximately 50~70 MHz can 
be assigned to IPs, although we cannot resolve three in- 
equivalent IPs precisely from the present ZF-NMR spec- 
tra. The Mafm at those IPs may be tentatively es- 
timated as Mafm (IPs) =0.24 0.28/^6 • The ZF-NMR 
spectra at approximately 120 ^ 150 MHz arise from 
the AFM-Mott insulator CaCu02(IL), which is contam- 
inated during synthesis process under high-pressure and 
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Fig. 19. (Color online) T dependence of K^^ at OP in nearly 
optimally doped Hgl2(n-l)?i series: (a) n=4:Hgl234(OPT) with 
rc=123 K,3'' (b) n=5:Hgl245(OPT)ttl with Tc=108 K,^^ and 
(c) n=8:Hgl278(OPT) with rc=103 K. The broken curves in (b) 
and (c) show the T dependence of at OP of Hgl234(OPT), 
which shows no AFM order. The Cu-NQR/ZFNMR spec- 
tra at Bcxt=0 and r=1.5 K for (d) n=4:Hgl234(OPT), (e) 
n=5:Hgl245(OPT)t)l,3'* and (f) n=8:Hgl278(OPT). The spec- 
tra indicated by the symbol (*) arise from the AFM-Mott in- 
sulator CaCu02(IL), which is contaminated during the high- 
pressure synthesis process. The inset at the bottom illustrates 
the physical properties of ri=8:Hgl278(OPT) with Tc=103 K. 
[cited from rcf.^^^] 



high-temperature conditions. Further extensive NMR 
measurements are desired for cuprates with more than 
six layers over a wide doping level. 

Nevertheless, note that T^—^-^i K for 
n=8:Hgl278(OPT) remains high, even though the 
OPs that are responsible for HTSC are separated by 
a thick AFM block consisting of six IPs. When noting 
that SC Cooper pairs can tunnel between OPs through 
AFM ordered IPs by virtue of the Josephson coupling 
as predicted by a theoretical study of multilayered 
cuprates, ^^■^ it would be expected that HTSC can 
maintain relatively high and constant values, as 
shown in Fig. 2(b). 

More importantly, when noting that the OP with 
p(OP)^ 0.14 is paramagnetic, we remark that the criti- 
cal hole density of n=8 compounds Pc(8) does not reach 
0.14 even when an interlayer magnetic coupling is suffi- 
ciently increased. This result is consistent with the fact 
that Mafm(p) at 1.5 K decreases monotonically as p 
increases and is extrapolated to zero in the range of 
Pc(Mafm=0) < 0.14, as discussed in § 4.2 (see Fig. 22). 
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Fig. 20. (Color online) Phase diagrams of AFM and HTSC for n=l:LSCO,3'*' n=2:0212F''S and YBCOe+rr/^^' ^^"^ n=3:0223F,'**' ^0 
ri=4:0234F,''l''"''''5 and n=5:M1245.3*' ""^^ Pc(n), at which an AFM order collapses, increases from pc{n) ~0.075 to 0.08 to 
0.10 as n increases from 3 to 4 to 5, respectively. The result on ?i=8 compound has revealed that Pc(8) does not reach 0.14 even 
when an interlayer magnetic coupling is enhanced, which is also deduced from the fact that Mafm decreases to zero in the range of 
Pc(Mafm=0) < 0.14 (see Fig. 22). [cited from refs.**^' ^2, 53] 



4. Discussion 

4-.1 Number n of Cu02 layer dependence of phase dia- 
gram of AFM and HTSC 
Figure 20 shows the phase diagrams of AFM 
and HTSC for n=l: LSCO,34.i24 „^2: 0212F49 and 
YBC06+:.,i25,i26 „=3:0223F,44'50 n=4:0234F,4i'44,45 

and n=5:M12A5.^^'^^''^^''^^'^^ The phase diagram on n- 
layered cuprates with n=3, 4, and 5 are characterized as 
follows: (i) The AFM metallic phase is robust up to an 
optimally hole-doped region: The quantum critical point 
(pc) at which an AFM order collapses increases from 
Pcin) ~0.075 to 0.08 to 0.10 as n increases from 3 to 4 to 

5, respectively, (ii) The uniform coexistence of AFM and 
HTSC is the universal phenomenon inherent in a single 
Cu02 plane in an underdoped region, (iii) The maximum 
Tc takes place at approximately p(T™*'^) ^ 0.16 irrespec- 
tive of n. 

The phase diagrams of n- layered cuprates with rt=3, 
4, and 5 differ significantly from the well-established 
phase diagrams of n=l:LSCO and n=2:YBC06+a:, where 



Pc{l) - 0.0234 p^^2) ~ 0.055,125.126 respectively, as 
shown in Fig. 20. In fact, from the present NMR mea- 
surement, the phase diagram of n=2:0212F does not re- 
veal an AFM order in the range of p > 0.083.^^ It is 
apparent that as n increases from 1 to 5, Pc{n) increases 
from 0.02 to 0.10. The mother compounds for HTSC are 
characterized by a large in-plane superexchange inter- 
action Jin ~1300 K among spins at nearest-neighbor Cu 
sites. ^"1° However, since no long-range AFM order occurs 
at a finite temperature for an isolated two-dimensional 
(2D) system, the interlayer magnetic coupling along the 
c-axis, which depends on n, plays a crucial role in sta- 
bilizing an AFM order. An effective interlayer magnetic 
coupling of n-layered cuprates is given as \/ JcJoutin), 
where Jc is the magnetic coupling between OPs through 
CRL and Jout is the magnetic coupling in a unit cell, 
as illustrated in Fig. 21(a). Here, Jc is independent of n, 
but Jout('^) increases with increasing n. In this context, it 
is the weak interlayer magnetic coupling that suppresses 
the static long-range AFM order in LSCO and YBCOg+x 
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Fig. 21. (Color online) Schematic illustrations of magnetic cou- 
plings in n-layered cuprates: (a) interlayer magnetic coupling 
along c-axis and (b) in-planc superexchange interaction Jj^ 
among spins at nearest-neighbor Cu sites in two-dimensional lat- 
tice of Cu02 plane. Here, Jc is the magnetic coupling between 
OPs through CRL, which is independent of n, and Jout(i) is the 
magnetic coupling between OPs through IPs, which increases 
with n. Since Jin is as large as 1300 K in undoped AFM-Mott 
insulators, an AFM order is stabilized when the interlayer 
magnetic coupling (\/ JcJout{n)) becomes stronger with increas- 
ing 71. 



at such small carrier densities. In the contrast, the re- 
sult for the n=8 compound has revealed that Pc(8) does 
not reach 0.14 even when the interlayer magnetic cou- 
pling becomes sufBciently large. It is likely that pc sat- 
urates at approximately 0.14 even in the strong limit of 
interlayer magnetic coupling expected for the n=oo com- 
pound, since Mafm in the ground state is extrapolated 
to zero at Pc(AfAFM=0) < 0.14, as discussed in the next 
section (see Fig. 22). 

4-2 Ground-state phase diagram of AFM and HTSC 

Since Tn depends on the strength of interlayer mag- 
netic coupling, i.e., the number of Cu02 layers, the tem- 
perature phase diagram against p depends on n, as shown 
in Fig. 20. Thus, we discuss here a ground-state phase 
diagram of a Cu02 plane by plotting the AFM mo- 
ment (AfAFM) and SC gap (Asc) against p, which also 
gives us an opportunity to compare the experimental out- 
comes with theoretical ones. 

Figure 22 shows the p dependence of the Mafm 
at T=1.5 K and Asc evaluated using the relation 
2Asc=8fcBTc.''3-46, 52-54 rj,^^ ]^^^^ ^^le CuOa plane 

decreases monotonically as p increases and is extrapo- 
lated to zero in the range of 0.12< Pc{Mafm—0) <0.14, 
as shown by the solid curve in Fig. 22. Here, the data 
of Mafms for nondoped and slightly doped Mott insu- 
lating states are cited from n=oo:Cao.85Sro.i5Cu02^''^ 
and n=2:YBC0a;,"^^^ respectively. Figures 20 and 22 re- 
veal that Pc(Mafm=0) ^0.14 is an intrinsic quantum 
critical hole density at the Cu02 plane, at which the 
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Fig. 23. (Color online) Plot of Tn vs Mafm- The 
solid curve shows Tn (X AflpM witti Tn=537 K at 
^AFM=0.51/xb in '^=c»:Cao.85Sro.i5Cu02.^''® The arrow points 
to Mafm ~ 0.18/tB- [cited from refs.''*''52,53] 



AFM moment totally disappears in the ground state. 
It is noteworthy that the maxima of Asc and Tc are 
at p(r™'*'^) 0.16 irrespective of n, which are close to 
Pc{Mafm=0) ^0.14, revealing an intimate relationship 
between the SC order parameter and the AFM ordered 
moment. The phase diagram of Mafm vs p presented 
here is totally consistent with the ground-state phase 
diagram at a single Cu02 plane, which has been the- 
oretically addressed thus far in terms of either the t-J 
model,"-!^' 19-22, 24-28, 31. 32 ^j^g Hubbard model in a 

strong correlation regime. Note that Pc(Mafm=0) ^ 
0.14 determined experimentally is also in good agreement 
with the theoretical calculation at T=0. 

We also point out that no static AFM order with a 
3D long-range character can be observed in the range of 
Pc{n) < p < Pc(Mafm=0) owing to the strong 2D fluctu- 
ations and weak interlayer magnetic coupling; hence, the 
AFM moments at Cu02 planes may fluctuate on a time 
scale faster than the nanosecond order in NMR measure- 
ment, which may be related to various anomalous mag- 
netic behaviors in underdoped regions, as addressed in § 
4.6. 

4-3 Superexchange interaction in hole-doped Cu02 
plane 

Figure 23 shows a plot of Tn vs A/afm to 
gain insight into the p dependence of the in- 
plane superexchange interaction Jin(p)- In this fig- 
ure, data are presented with respect to n=3:0223F,^'* 
n=4:0234F,45 n=^:M\2Ah,'^^'^^ along with the data 
of n=oo:Cao.85Sro.i5Cu02 with p—Q}'^^ 

On the basis of the mean-field approximation of local- 
ized spins, we assume that oc M^pj^, and Jout("- — 
oo) and Jin(p = 0) for Cao.85Sro.i5Cu02 stay constant 
regardless of p. The solid curve in Fig. 23 shows a 
graph of Tjv oc Mj^-py^ plotted using TAr(oo)=537 K 
at Mafm=0.51 /iB for Cao.85Sro.i5Cu02. First, we no- 
tice that, as shown by an arrow in Fig. 23, T^in) 
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Fig. 22. (Color online) Ground-state phase diagram: p-dependence of Mafm at T=1.5 K and SC gap (Age) for n=2:0212F,*® 
n=3:0223F,'"' n=4:0234F,''5 and n=5:M1245.''3' Data of Mafms for nondoped and slightly doped Mott insulating states are 
cited from n=oo:Cao.85Sro.i5Cu02^'^® and n=2:YBC06+a:/^^ respectively. The AFM moment at the Cu02 plane totally disappears 
in the ground state when p= 0.12~0.14 (= Pc(Mafm=0)), which is extrapolated from the p dependence of AfAFM(solid curve). Age 
shows a maximum just outside of Pc(Mafm=0)~0.14 irrespective of n. Here, the Age is estimated from the Tc values with the relation 
2Asc=8/cbT'c. [cited from refs.^^. 52-54] 



increases from n—A to 5 as n increases, even though 
-^AFM ~ 0.18 /xb remains constant for the n=4 and 
5 compounds, which is attributed to the increase in 
-'out('^), namely, Jout(4) < Jout(5). The effective in- 
terlayer couphng of the n—A and 5 compounds, given 
by a/ JcJoutin), is always smaller than Jout(oo) in 
Cao.85Sro.i5Cu02: However, as shown by the solid curve 
in Fig. 23, Tiy{n)s for the n=4 and 5 compounds, which 
are given by Tjyin) ~ Mlp^{p)[Jin{p)^y JcJout{n)]^^'^ , 
are larger than Tn{oo) for n=oo:Cao.85Sro.i5Cu02, 
which is given by TAr(oo) Af^pM(p)[Jin(0) Jout(oo)]i/2 
with T7v=537 K at MAFM=0.51/iB and p=0. Thus, 
we also obtain an unexpected relation, i.e., Jin(p) > 
</in(0) ~ 1300 K. The two experimental relationships 
— the plot of Mafm vs p shown in Fig. 22 and the 
plot of Tn vs Mafm shown in Fig. 23 — suggest that 
the AFM ground state in the homogeneously hole-doped 
Cu02 planes is determined by p, \/ JcJout{n), and Jin(p) 
that is larger than Jin(O) ^ 1300 K. It is surprising 
that Jin(p) becomes stronger for the doped Cu02 planes 
with AFM order than for the undoped AFM-Mott in- 
sulators. Mean-field theories of HTSC used to consider 
the superexchange interaction Jin as the source of an in- 
stantaneous attraction that led to pairing in a c?-wave 
state. The present outcomes may experimentally sup- 
port such a scenario as far as the underdoped region is 
concerned, where AFM and HTSC uniformly coexist at 
a Cu02 plane. Further theoretical study is desired to ad- 
dress whether Jin (p) becomes larger than that for AFM- 
Mott insulators even though mobile holes are doped. 



4-4 Pseudo gap behavior in underdoped region 

The pseudogap behavior emerging above Tc is the un- 
derlying issue in cuprate superconductors. NMR and 
ARPES studies have observed pseudogap behaviors in 
spin and charge excitations below r*pj„ and T*f^^^.g^, re- 
spectively. Figure 24 shows the p dependence of T* for 
n=5:M1245, along with T* for n=l:Bi220l"^' 
Hgl201,i2^ n=2:Bi2212,^«'93,i29 n=3:0223F5" and 
M1223.^^'^^° As shown in Fig. 24, the p dependences 
of T*pj„ and Tcharge rescmblc each other for the n=l 
and n—2 compounds, as argued in the literature. "'^^^ Note 
that, as far as no AFM order does not takes place, the 
p dependence of T/pj„ in the n=5 compounds resembles 
those for the n—1, 2, and 3 compounds, suggesting that 
the magnetic and charge excitations are suppressed be- 
low the same T, that is, T*^^^ coincides with T*^^^^^^. 
Despite the stronger magnetic interlayer coupling in the 
n=5 compounds than in the other multilayered cuprates, 
the p-dependence of T*^^^ resembles those of n=l and 
n—2 compounds in the optimally doped regime: it is ex- 
pected that T* is determined by in-plane magnetic and 
charge correlations. "'^'^"'^ From a theoretical point of view, 
the t-J model has explained T* as an onset tempera- 
ture for a singlet pairing of spinons in the Cu02 plane, '^^ 
which would be insensitive to the presence of magnetic 
interlayer coupling. 

Here, we highlight that the spin-gap behavior disap- 
pears when the AFM order takes place in p < Pc{ri); 
l/TiT continues to increase toward Tn upon cooling for 
the n=3 and n=b compounds that exhibit the AFM or- 
der at low temperatures, as shown in Figs. 11 and 18, 
respectively. In underdoped regions where the ground 
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Fig. 24. (Color online) p dependences of T* for n=5:M1245, 
along with the T* for n=l:Bi220l"*' and Hgl201,i2* 
n=2:Bi2212 93, 129 n=3:0223F50 and M1223.'^^' "° Those 
data were obtained by the measurements of ARPES^'^^' and 
NMR.^"'*^^'^*' 93. 118. 128, 130 rpj^g broken curve represents the p 
dependence of Tc, Tc=r™''''[l-82.6(p-0.16)2].90, 132 ^ 
AFM order takes place, the p dependences of T*^^^ and T'*^^^,^^ 
resemble irrespective of n, suggesting that T* is determined by 
in-plane magnetic and charge correlations, [cited from ref.^"] 



state is characterized by the coexistence of HTSC and 
AFM, low-lying magnetic excitations develop upon cool- 
ing toward the AFM order without a gap opening. In 
this context, unless an AFM order occurs in the under- 
doped region, it is natural for spin-singlet formation to 
develop above leading to a pseudogap behavior in 
magnetic excitations. Recently, Yamase et al. have the- 
oretically pointed out that a spin gap is strongly sup- 
pressed near a tetracritical point for phases of the AFM, 
HTSC, AFM-hHTSC, and normal states.^^ As a resuh, 
it is likely that the pseudogap, which emerges around 
the antinode region at the wave vectors (0,7r) and (7r,0), 
evolves to a real AFM gap in the AFM-HTSC mixed 
phase. One underlying issue is to experimentally address 
whether or not a pseudogap survives even in the para- 
magnetic and normal states in the AFM-HTSC mixed 
region. 

^.5 Momentum dependences of SC gap and pseudogap 
in HTSC 

As observed in ARPES experiments on multilayered 
cuprates,^^^' the momentum dependence of the gap 
magnitude in HTSC is schematically drawn in Fig. 25. 
Here, the gap for overdoped HTSC is almost a simple d 
wave, AqI cos(fc2;a) — cos(A:j,a)|/2, as shown by a straight 
line. On the other hand, the gap for the underdoped one 
deviates from the simple d-wave around the antinodes 
~ (0,7r) and (tt, 0). The gap size is characterized by two 
parameters, i.e., Aq around the node and A* in the antin- 
odal region, where Aq and A* are defined by the lin- 
ear extrapolation of the gap magnitude to the antinode 
[Ao| cos(fca;a) — cos(fcj,a)|/2 = Aq], as shown in the figure. 



underdope 
' overdope 




1.0 



Fig. 25. (Color online) Illustration of momentum dependences 
of energy gap Aq and pseudogap A* .1'^'' [cited from ref.^^] 



The deviation of the gap anisotropy from the simple d- 
wave is known to be prominent in underdoped cuprates, 
which is called the two-gap behavior .^^"^ 

The ARPES on n=4:0234F(jl4) with Tc=55 K showed 
two Fermi sheets of IP and OP and revealed that their 
momentum dependences of the gap magnitude exhibit a 
two-gap behavior at IP and OP at 20 K.^^^ When incor- 
porating that the AFM order sets in at Tn^SO K well 
above rc=55 K and that the AFM unit cell doubles a 
crystal one. A* around the antinodes ^ (0, tt) and (tt, 0) 
is assigned to an AFM gap. By contrast. A*, which used 
to be observed even in the absence of an AFM order in 
the underdoped regime for the n=\ and 2 compounds, 
may be ascribed to a formation of a spin gap due to the 
development of a singlet resonating valence bond (RVB) 



state 



31 



^.6 Comments on underlying issues in underdoped 
cuprates 

The experimental results in Figs. 20 and 22 sug- 
gest that an AFM moment exists at Cu sites up 
to Pc(^AFM=0) < 0.14 even in the n=l:LSCO and 
n=2:YBC06+2: compounds, but it cannot be observed 
as a static AFM order in the range of Pc{n) < p < 
Pc{Mafm=0) owing to strong 2D fluctuations since the 
interlayer magnetic coupling is too weak to stabilize a 
completely static AFM order with a 3D long-range char- 
acter. This instability of the ground state on the elec- 
tronic states was corroborated by the facts that the ap- 
plication of a high magnetic field helps to bring about 
an AFM order with Mafm ~ 0.1 /LtB in the vortex state 
of LSCO with p=0.li35 and ¥6006.5^^^ and that the 
spin-glass phase existing in LSCO survives up to a hole 
density close to Pc{Mafm=0) ~ 0.14.^^"^ In addition, the 
static AFM order was also observed when the charge- 
stripe order occurs at x ~l/8 of LSCO.^"^^ The neutron 
scattering experiment suggested that an unusual AFM 
order with Mafm ~ O.I/ib takes place in the under- 
doped YBCOg.s with Tc—60 K, although it fluctuates on 
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a nanosecond time scale. The phase diagrams of LSCO 
and YBCOe+x have been beheved as the prototypes thus 
far; however, we note that the underlying issues in the 
underdoped region mentioned above should be reconsid- 
ered on the basis of the following facts: (i) an AFM mo- 
ment intrinsically exists up to Pc{Mafm=0) < 0.14 in 
the ground state, (ii) the intcrlaycr magnetic coupling 
in LSCO and YBCOe+x is not sufficient to stabilize a 
completely static AFM order and (iii) the chemical sub- 
stitution for doping introduces inevitable disorders into 
the Cu02 plane. 

5. Concluding Remarks 

Site-selective NMR studies on multilayered cuprates 
have unraveled the intrinsic phase diagram of homoge- 
neously hole-doped Cu02 plane as follows: 

(i) The AFM metallic phase is robust and uniformly 
coexists with the HTSC phase up to the quantum 
critical hole density {pc{n)), at which the AFM order 
collapses. 

(ii) Pc{n) increases from 0.075 to 0.08 to 0.10 as the in- 
terlayer magnetic coupling becomes stronger in in- 
creasing from n=3 to 4 to 5, respectively. 

(iii) The uniform coexistence of AFM and HTSC at a 
single Cu02 plane is the universal phenomenon in 
the underdoped region for p < Pc{n). 

(iv) The maxima of Asc and Tc take place irrespective 
of n at piT^'''') ~ 0.16 just outside Pc{Mafm=0) - 
0.14, at which the AFM moment inherent in the 
Cu02 plane disappears in the ground state. 

(v) The pseudogap in magnetic excitations collapses in 
the AFM state for p < Pc{n), where low- lying mag- 
netic excitations develop upon cooling toward the 
AFM order. 

(vi) The ground-state phase diagram of AFM and HTSC 
(see Fig. 22) is in good agreement with the ground- 
state phase diagrams in terms of either the t-J 
model"-!^' 19-22, 24-28, 31, 32 ^^iB Hubbard model in 
a strong correlation regime. 3" 

(vii) The in- plane superexchange interaction Jiaip) 
for the n=4 and 5 compounds is larger than 
Ji„(0) - 1300 K for the infinite-layered AFM-Mott 
insulator Cao.85Sro.i5Cu02. 

In particular, we emphasize that the uniformly mixed 
phase of AFM and HTSC for p < Pc{n) and the emer- 
gence of d-wave SC with the maximum of Tc just outside 
Pc{M AFM=0) can be accounted for by the Mott physics 
based on the t-J model. In the physics behind high-Tc 
phenomena, there is a very strong Coulomb repulsive in- 
teraction ?7 (> 6 eV), which prohibits the double occu- 
pancy of an up-spin electron and a down-spin electron on 
the same site. When noting that U is almost unchanged 
with doping, being almost the same as that in AFM-Mott 
insulators, a large Jjn attracts electrons of opposite spins 
to be in neighboring sites, raising the Tc of cuprates up 
to as high as 160 K;^^, i3,3i ^^ere are no bosonic glues. 



Acknowledgement 

These works have been carried out in collaboration 
with M. Abe, Y. Yamaguchi, T. Sakaguchi, K. Itohara, 
S.-i. Tabata, S. Iwai, K. Matoba, Y. Araki, H. Kotegawa, 
Y. Tokunaga, K. Magishi, K. Ishida, G.-q. Zheng, and 
K. Asayama for NMR studies of multilayered cuprates. 
The samples were provided by P. M. Shirage, H. Kito, 
Y. Kodama, Y. Tanaka, H. Ihara (AIST), K. Tokiwa, 
and T. Watanabe (Tokyo University of Science). We 
thank M. Mori, T. Tohyama, S. Maekawa, H. Yamase, 
T. K. Lee, M. Ogata, and H. Fukuyama for their valu- 
able discussions and comments. These works were sup- 
ported by a Grant-in-Aid for Specially Promoted Re- 
search (20001004) and by the Global COE Program 
(Core Research and Engineering of Advanced Materials- 
Interdisciplinary Education Center for Materials Science) 
from the Ministry of Education, Culture, Sports, Science 
and Technology (MEXT), Japan. 



6 



J. G. Bcdnorz and K. A. MuUer; Z. Phys. B 64 (1986) 189. 
A. Schilling, M. Cantoni, J. D. Guo, and H. R. Ott: Nature 

(London) 363 (1993) 56. 

L. Gao, Y. Y. Xuc, F. Chen, Q. Xiong, R. L. Mcng, D. Ramirez, 
and C. W. Chu: Phys. Rev. B 50 (1994) 4260. 
C. W. Chu, L. Gao, Z. J. Huang, R. L. Meng, and Y. Y. Xue: 
Nature (London) 365 (1993) 323. 

H. Ihara, M. Hirabayaslii, H. Tanino, K. Tokiwa, H. Ozawa, Y. 
Akahama, and H. Kawamura: Jpn. J. Appl. Phys. 32 (1993) 

L1732. 

J. Bardccn, L. N. Cooper, and J. R. Schriefler: Phys. Rev. 108 
(1957) 1175. 

S. M. Hayden, G. Aepph, R. Osborn, A. D. Taylor, T. G. Pcr- 
rlng, S.-W. Cheong, and Z. Fisk: Phys. Rev. Lett. 67 (1991) 
3622. 

K. B. Lyons, P. A. Fleury, L. F. Schneemeyer, and J. V. 
Waszczak: Phys. Rev. Lett. 60 (1988) 732. 
S. M. Hayden, G. Acppli, T. G. Pcrring, H. A. Mook, and F. 
Dogan: Phys. Rev. B 54 (1996) R6905. 

Y. Tokura, S. Koshihara, T. Arima, H. Takagi, S. Ishibashi, 
T. Ido, and S. Uchida: Phys. Rev. B 41 (1990) 11657. 
P. W. Anderson: Science 235 (1987) 1196. 
P. W. Anderson; The Theory of Superconductivity in the High- 
Tc Cuprate Superconductors (Princeton Univ. Press, Prince- 
ton, 1997). 

P. W. Anderson, P. A. Lee, M. Randeria, T. M. Rice, N. Trivedi, 
and F. C. Zhang: J. Phys. Condcns. Matter 16 (2004) R755. 
G. J. Chen, R. ,Ioynt, F. C. Zhang, and C. Gros: Phys. Rev. 
B 42 (1990) 2662. 

T. Giamarchi and C. Lhuillier; Phys. Rev. B 43 (1991) 12943. 
M. Inaba, H. Matsukawa, M. Saitoh, and H. Fukuyama: Phys- 
ica C 257 (1996) 299. 

T. K. Lee and C.T. Shih: Phys. Rev. B 55 (1997) 5983. 
S. C. Zhang: Science 275 (1997) 1089. 
A. Himeda and M. Ogata: Phys. Rev. B 60 (1999) R9935. 

G. Kotliar, S. Savrasov, G. Palsson, and G. Biroli: Phys. Rev. 
Lett. 87 (2001) 186401. 

A. Paramckanti, M. Randeria, and N. Trivedi: Phys. Rev. Lett. 
87 (2001) 217002. 

T. K. Lee, C. M. Ho, and N. Nagaosa: Phys. Rev. Lett. 90 

(2003) 067001. 

E. Demler, W. Hanke, and S. C. Zhang: Rev. Mod. Phys. 76 

(2004) 909. 

H. Yamase and H. Kohno: Phys. Rev. B 69 (2004) 104526. 
H. Yamase, M. Yoneya, and K. Kuboki: Phys. Rev. B 84 (2011) 
014508. 

26) A. Paramckanti, M. Randeria, and N. Trivedi: Phys. Rev. B 
70 (2004) 054504. 



13 
14 

15 
16^ 

17 

18; 

19^ 
20 

21 

22 

23 

24 
25' 



J. 



Phys. Soc. Jpn. 



Full Paper 



Author Name 19 



C. T. Shih, T. K. Lcc, R. Edcr, C. Y. Mou, and Y. C. Chen: 60 
Phys. Rev. Lett. 92 (2004) 227002. 

C. T. Shih, J. J. Wu, Y. C. Chen, C. Y. Mou, C. P. Chou, R. 61 
Eder, and T. K. Lee: Low Temp. Phys. 31 (2005) 757. 

D. Senechal, P.-L. Lavertu, M.-A. Marois, and A.-M. S. Trem- 62 
blay: Phys. Rev. Lett. 94 (2005) 156404. 

M. Capone and G. KotUar: Phys. Rev. B 74 (2006) 054513. 63 
M. Ogata and H. Pukuyama: Rep. Prog. Phys. 71 (2008) 

036501. 

S. Pathak, V. B. Shcnoy, M. Randeria, and N. Trivedi: Phys. 64 
Rev. Lett. 102 (2009) 027002. 

T. Moriya and K. Ueda: Rep. Prog. Phys. 66 (2003) 1299. 65 
B. Keimer, N. Belk, R. J. Birgeneau, A. Cassanho, C. Y. Chen, 
M. Greven, M. A. Kastner, A. Aharony, Y. Endoh, R. W. Erwin, 
and G. Shirane: Phys. Rev. B 46 (1992) 14034. 66 
Y. Tokunaga, H. Kotegawa, K. Ishida, G.-q. Zheng, Y. Kitaoka, 
K. Tokiwa, A. lyo, and H. Ihara: J. Low Temp. Phys. 117 (1999) 67 
473. 

Y. Tokunaga, K. Ishida, Y. Kitaoka, K. Asayama, K. Tokiwa, 68 
A. lyo, and H. Ihara; Phys. Rev. B 61 (2000) 9707. 
H. Kotegawa, Y. Tokunaga, K. Isliida, G.-q. Zheng, Y. Kitaoka, 69 
K. Asayama, H. Kito, A. lyo, H. Ihara, K. Tanaka, K. Tokiwa, 
and T. Watanabe: Phys. Rev. B 64 (2001) 064515. 70 
H. Kotegawa, Y. Tokunaga, Y. Araki, G.-q. Zheng, Y. Kitaoka, 
K. Tokiwa, K. Ito, T. Watanabe, A. lyo, Y. Tanaka, and H. 71 
Ihara; Phys. Rev. B 69 (2004) 014501. 

H. Mukuda, M. Abe, Y. Araki, Y. Kitaoka, K. Tokiwa, T. 
Watanabe, A. lyo, H. Kito, and Y. Tanaka: Phys. Rev. Lett 72 

96 (2006) 087001. 

H. Mukuda, M. Abe, S. Shimizu, Y. Kitaoka, A. lyo, Y. Ko- 
dama, H. Kito, Y. Tajiaka, K. Tokiwa, and T. Watanabe: J. 73 
Phys. Soc. Jpn. 75 (2006) 123702. 74 
S. Shimizu, H. Mukuda, Y. Kitaoka, A. lyo, Y. Tanaka, Y. 
Kodama, K. Tokiwa, and T. Watanabe; Phys. Rev. Lett. 98 75 
(2007) 257002. 76; 
H. Mukuda, M. Abe, Y. Kitaoka, H. Kotegawa, K. Tokiwa, 
T. Watanabe, A. lyo, H. Kito, Y. Tanaka, and Y. Kodama: 77 
Physica C 460-462 (2007) 36. 

H. Mukuda, Y. Yamaguchi, S. Shimizu, Y. Kitaoka, P. M. Shi- 
rage, and A. lyo: J. Phys. Soc. Jpn. 77 (2008) 124706. 78 
S. Shimizu, T. Sakaguchi, H. Mukuda, Y. Kitaoka, P. M. Shi- 
rage, Y. Kodama, and A. lyo: Phys. Rev. B 79 (2009) 064505. 79 
S. Shimizu, H. Mukuda, Y. Kitaoka, H. Kito, Y. Kodama, P. 
M. Shiragc, and A. lyo: J. Phys. Soc. Jpn. 78 (2009) 064705. 80 
H. Mukuda, S. Shimizu, S.-i. Tabata, K. Itohara, Y. Kitaoka, 
H. Kito, P. M. Shiragc, and A. lyo; Physica C 470 (2010) S7. 81 
K. Itohara, S. Shimizu, H. Mukuda, Y. Kitaoka, P. M. Shiragc, 
H. Kito, and A. lyo: Physica C 470 (2010) S140. 82 
S. Shimizu, S.-i. Tabata, H. Mukuda, Y. Kitaoka, P. M. Shirage, 
H. Kito, and A. lyo: J. Phys. Soc. Jpn. 80 (2011) 043706. 83 
S. Shimizu, S. Iwai, S.-i Tabata, H. Mukuda, Y. Kitaoka, P. M. 
Shirage, H. Kito, and A. lyo: Phys. Rev. B 83 (2011) 144523. 84 
S. Shimizu, S.-i. Tabata, H. Mukuda, Y. Kitaoka, P. M. Shirage, 85 
H. Kito, and A. lyo: Phys. Rev. B 83 (2011) 214514. 
S.-i. Tabata, S. Shimizu, H. Mukuda, Y. Kitaoka, P. M. Shirage, 
and A. lyo: unpublished. 86 
Y. Kitaoka, H. Mukuda, S. Shimizu, S.-i. Tabata, P. M. Shirage, 87 
and A. lyo; J. Phys. Chem. Solids 72 (2011) 486. 
Y. Kitaoka, S. Shimizu, H. Mukuda, S.-i. Tabata, P. M. Shirage, 
and A. lyo: J. Phys. Condens. Matter (2011) in press. 
S. Shimizu, S.-i. Tabata, S. Iwai, H. Mukuda, Y. Kitaoka, P. 
M. Shirage, H. Kito, and A. lyo: Phys. Rev. B in press. 89 
K. Tokiwa, A. lyo, T. Tsukamoto, and H. Ihara: Czech. J. Phys. 
46 (1996) 1491. 90 
K. Tokiwa, S. Ito, H. Okumoto, S. Mikusu, A. lyo, Y. Tanaka, 
and T. Watanabe; J. Low. Temp. Pliys. 131 (2003) 637. 91 
A. lyo, Y. Tanaka, H. Kito, Y. Kodama, P. M. Shirage, D. 
D. Shivagan, H. Matsuhata, K. Tokiwa, and T. Watanabe: J. 92 
Phys. Soc. Jpn. 76 (2007) 094711. 

A. lyo, Y. Tanaka, M. Tokumoto, and H. Ihara: Physica C 366 
(2001) 43. 93 
A. lyo, M. Hirai, K. Tokiwa, T. Watanabe, and Y. TanaJja: 
Physica C 392-396 (2003) 140. 



A. lyo, Y. Tanaka, Y. Kodama, H. Kito, K. Tokiwa, and T. 
Watanabe: Physica C 445 (2006) 17. 

P. M. Shirage, D. D. Shivagan, Y. Tanaka, Y. Kodama, H. Kito, 
and A. lyo: Appl. Phys. Lett. 92 (2008) 222501. 
M. Hirai, Y. Kodama, H. Kito, Y. Tanaka, K. Tokiwa, T. 
Watanabe, and A. lyo: J. Phys. Soc. Jpn. 76 (2007) 054701. 

B. A. Scott, E. Y. Suard, C. C. Tsuei, D. B. Mitzi, T. R. 
McGuire, B.-H. Chen, and D. Walker: Physica C 230 (1994) 
239. 

E. V. Antipov, A. M. Abakumov, and S. N. Putilin: Supercond. 
Sci. Technol. 15 (2002) R31. 

M.-H. Julien, P. Carrctta, M. Horvatic, C. Berthier, Y. 
Berthier, P. Segransan, A. Carrington, and D. Colson: Phys. 
Rev. Lett. 76 (1996) 4238. 

K. Magishi, Y. Kitaoka, G.-q. Zheng, K. Asayama, K. Tokiwa, 

A. lyo, and H. Ihara: J. Phys. Soc. Jpn. 64 (1995) 4561. 

K. Magishi, Y. Kitaoka, G.-q. Zheng, K. Asayama, K. Tokiwa, 

A. lyo, and H. Ihara: Phys. Rev. B 53 (1996) R8906. 

S. Ohsugi, Y. Kitaoka, K. Ishida, G.-q. Zheng, and K. Asayama: 

J. Phys. Soc. Jpn. 63 (1994) 700. 

Y. Itoh, T. Machi, S. Adachi, A. Fukuoka, K. Tanabe, and H. 

Yasuoka: J. Phys. Soc. Jpn. 67 (1998) 312. 

K. Ishida, Y. Kitaoka, K. Asayama, K. KadowaJci, and T. 

Mochiku: J. Phys. Soc. Jpn. 63 (1994) 1104. 

S. E. Barrett, D. J. Durand, C. H. Pennington, C. P. Slichter, 

T. A. Fricdmann, J. P. Rice, and D. M. Ginsberg: Phys. Rev. 

B 41 (1990) 6283. 

H. Eisaki, N. Kancko, D. L. Feng, A. Damascclli, P. K. Mang, K. 
M. Shen, Z.-X. Shen, and M. Greven: Phys. Rev. B 69 (2004) 
064512. 

F. Mila and T. M. Rice: Phys. Rev. B. 40 (1989) 11382. 

A. J. Millis, H. Monien, and D. Pines: Phys. Rev. B 42 (1990) 
167. 

T. Imai: J. Phys. Soc. Jpn. 59 (1990) 2508. 

H. Monien, D. Pines, and M. Takigawa: Phys. Rev. B 43 (1991) 

258. 

M. Takigawa, A. P. Reyes, P. C. Hamnicl, J. D. Thompson, R. 
H. Heffner, Z. Fisk, and K. C. Ott: Piiys. Rev. B 43 (1991) 
247. 

R. E. Walstedt, R. F. Bell, and D. B. Mitze: Phys. Rev. B. 44 
(1991) 7760. 

H. Zimmermann, M. Mali, M. Bankay, and D. Brinkmann: 
Physica C 185-189 (1991) 1145. 

M. W. Shafcr, T. Penney, and B. L. Olson: Phys. Rev. B 36 

(1987) 4047. 

K. Kishino, J.-i. Shimoyama, T. Hascgawa, K. Kitaaawa, and 
K. Fucki: Jpn. J. Appl. Phys. 26 (1987) L1228. 
Y. Tokura, J. B. Torrance, T. C. Huang, and A. I. Nazzal: Phys. 
Rev. B 38 (1988) 7156. 

J. B. Torrance, Y. Tokura, A. I. Nazzal, A. Bezinge, T. C. 
Huang, and S. S. P. Parkin: Phys. Rev. Lett. 61 (1988) 1127. 

I. D. Brown: J. Sohd State Chem. 82 (1989) 122. 

R. J. Cava, A. W. Hewat, E. A. Hewat, B. Batlogg, M. Marezio, 
K. M. Rabe, J. J. Krajewski, W. F. Peck Jr., and L. W. Rupp 
Jr.: Physica C 165 (1990) 419. 
J. L. Tallon: Physica C 168 (1990) 85. 

A. A. Kordyuk, S. V. Borisenko, M. S. Golden, S. Legner, K. 
A. Nenkov, M. Knupfer, J. Fink, H. Berger, L. Forro, and R. 
Follath: Phys. Rev. B 66 (2002) 014502. 
S. D. OberteUi, J. R. Cooper, and J. L. Tallon: Phys. Rev. B 
46 (1992) 14928. 

J. L. Tallon, C. Bernhard, H. Shaked, R. L. Hitterman, and J. 

D. Jorgensen: Phys. Rev. B 51 (1995) R12911. 

M. R. Presland, J. L. Tallon, R. G. Buckley, R. S. Liu, and N. 

E. Flower: Physica C 176 (1991) 95. 

K. Fujiwara, Y. Kitaoka, K. Asayama, Y. Shimakawa, T. Man- 

ako, and Y. Kubo: J. Phys. Soc. Jpn. 59 (1990) 3459. 

R. E. Walstedt, W. W. Warren, Jr., R. F. Bell, R. J. Cava, 

G. P. Espinosa, L. F. Schneemeyer, and J. V. Waszczak: Phys. 
Rev. B 41 (1990) 9574. 

K. Ishida, K. Yoshida, T. Mito, Y. Tokunaga, Y. Kitaoka, K. 
Asayama, Y. Nakayama, J. Shimoyama, and K. Kishio: Phys. 
Rev. B 58 (1998) R5960. 



20 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



94) J. G. Storey, J. L. Tallon, and G. V. M. Williams: Phys. Rev. 
B 76 (2007) 174522. 

95) K. Magishi, Y. Kltaoka, G.-q. Zheng, K. Asayama, T. Kondo, 
Y. ShimaJjawa, T. Manako, and Y. Kubo: Phys. Rev. B 54 
(1996) 10131. 

96) G.-q. Zheng, Y. Kitaoka, K. Ishida, and K. Asayama: J. Phys. 
Soc. Jpn. 64 (1995) 2524. 

97) J. Haase, O. P. Sushkov, P. Horsch, and G. V. M. WiUiams: 
Phys. Rev. B 69 (2004) 094504. 

98) Y. Kitaoka, K. Pujiwara, K. Ishida, K. Asayama, Y. Shi- 
makawa, T. Manako, and Y. Kubo: Physica C 176 (1991) 
107. 

99) Y. Tanaka, A. lyo, N. Shirakawa, M. Ariyama, M. Tokumoto, 
S. I. Ikeda, and H. Ihara: Physica C 357-360 (2001) 222. 

100) A. Crisan, Y. Tanaka, A. lyo, L. Cosereanu, K. Tokiwa, and 
T. Watanabe: Phys. Rev. B 74 (2006) 184517. 

101) Y. Tanaka: Phys. Rev. Lett. 88 (2001) 017002. 

102) T. Tsuda, T. Olmo, and H. Yasuoka: J. Phys. Soc. Jpn. 61 
(1992) 2109. 

103) H. Yasuoka, T. Shimizu, Y. Ueda, and K. Kosugc: J. Phys. 
Soc. Jpn. 57 (1988) 2659. 

104) D. Vaknin, S. K. Sinha, D. E. Moncton, D. C. Johnston, J. M. 
Newsam, C. R. Saiinya, and H. E. King Jr.: Phys. Rev. Lett. 
58 (1987) 2802. 

105) D. Vaknin, E. Caignol, P. K. Davies, and J. E. Fischer: Phys. 
Rev. B 39 (1989) 9122. 

106) In the previous paper, the peak at ~160 MHz for 
Hgl245(UD)j)l was incorrectly assigned to IPs, which 
stemmed from some impurity phase. We reassign the spectra 
at 30~60 MHz to IPs of Hgl245(UD)ttl, which was detected 
under an appropriate pulse condition owing to the short T2.^^ 

107) K. Tokiwa, H. Okumoto, T. Imamura, S. Mikusu, K. Yuasa, W. 
Higemoto, K. Nishiyama, A. lyo, Y. Tanaka, and T. Watanabe: 
Int. J. Mod. Phys. B 17 (2003) 3540. 

108) K. Tokiwa, S. Mikusu, W. Higemoto, K. Nishiyama, A. lyo, Y. 
Tanaka, H. Kotcgawa, H. Mukuda, Y. Kitaoka, and T. Watan- 
abe: Physica C 460-462 (2007) 892. 

109) K. Nakayama and T. Moriya: J. Phys. Soc. Jpn. 56 (1987) 
2918. 

110) We deduce that SfnJF) in 0245F is dominated by a 
transferred hyperfine field from the nearest-neighboring 
Mafm(OP) not by Mafm(IP), since this Bf^^.{F) is observed 
along the c-axis owing to the possible out-of-plane canting of 
the Cu-3d spin direction at an antisymmetric Cu position of 
OP.« 

111) F. Borsa, P. Carrctta, J. H. Cho, F. C. Chou, Q. Hu, D. C. 
Jhonston, A. Lascialfari, D. R. Torgeson, R. J. Gooding, N. 
M. Salem, and K. J. E. Vos: Phys. Rev. B 52 (1995) 7334. 

112) H. Yasuoka, T. Imai, and T. Shimizu: in Strong Correlation 
and Superconductivity, ed. H. Pukuyama, S. Maekawa, and A. 
P. Malozemoff (Springer, Berlin, 1989) Vol. 89, p. 254. 

113) J. Rossat-Mignod, L. P. Regnault, C. Vettier, P. Burlet, J. Y. 
Henry, and G. Lapertot: Physica B 169 (1991) 58. 

114) A. G. Loeser, Z.-X. Shen, D. S. Dessau, D. S. Marshall, C. H. 
Park, P. Fournicr, and A. Kapitulnik: Science 273 (1996) 325. 

115) H. Ding, T. Yokoya, J. C. Campuzano, T. Takahashi, M. Ran- 
deria, M. R. Norman, T. Mochiku, K. Kadowaki, and J. Gi- 
apintzakis: Nature 382 (1996) 51. 

116) J. M. Harris, Z.-X. Shen, P. J. White, D. S. Marshall, M. C. 



Schabel, J. N. Eckstein, and I. Bozovic: Phys. Rev. B 54 (1996) 
R15665. 

117) R. E. Walstedt: The NMR Probe of High-Tc Materials 
(Springer, BerUn, 2008). 

118) G.-q. Zheng, P. L. Kuhns, A. P. Reyes, B. Liang, and C. T. 
Lin: Phys. Rev. Lett. 94 (2005) 047006. 

119) Y. Chen, A. lyo, W. Yang, X. Zhou, D. Lu, H. Eisaki, T. P. 
Devcrcaux, Z. Hussain, and Z.-X. Shen: Phys. Rev. Lett. 97 
(2006) 236401. 

120) W. Xic, O. Jcpscn, O.K. Andersen, Y. Chen, and Z.-X. Shen: 
Phys. Rev. Lett. 98 (2007) 047001. 

121) H. Mukuda, Y. Yamaguchi, S. Shimizu, Y. Kitaoka, P. M. 
Shirage, and A. lyo: unpublished. 

122) N. Kobayashi, Z. Hiroi, and M. Takano: J. Sol. State Chem. 

132 (1997) 274. 

123) M. Mori and S. Maekawa: Phys. Rev. Lett. 94 (2005) 137003. 

124) M.-H. Julien: Physica B 329-333 (2003) 693. 

125) S. Sanna, G. AUodi, G. Concas, A. D. Hillier, and R. De Renzi: 
Phys. Rev. Lett. 93 (2004) 207001. 

126) F. Concri, S. Sanna, K. Zliong, J. Lord, and R. De Renzi: 
Phys. Rev. B 81 (2010) 104507. 

127) T. Kondo, R. Khasanov, T. Takcuchi, J. Schmalian, and A. 
Kaminski: Nature 457 (2009) 296. 

128) Y. Itoh, T. Machi, A. Fukuoka, K. Tanabe, and H. Yasuoka: 
J. Phys. Soc. Jpn. 65 (1996) 3751. 

129) J. C. Campuzano, H. Ding, M. R. Norman, H. M. Eretwell, 
M. Randeria, A. Kaminski, J. Mesot, T. Takeuchi, T. Sato, 
T. Yokoya, T. Takahashi, T. Mochiku, K. Kadowaki, P. Gup- 
tasarma, D. G. Hinks, and Z. Konstantinovic: Phys. Rev. Lett. 
83 (1999) 3709. 

130) H. Kotcgawa, Y. Tokunaga, K. Ishida, G.-q. Zheng, Y. Ki- 
taoka, A. lyo, Y. Tanaka, and H. Ihara: Phys. Rev. B 65 (2002) 
184504. 

131) T. Yoshida, M. Hashimoto, S. Ideta, A. Fujimori, K. Tanaka, 
N. Mannella, Z. Hussain, Z.-X. Shen, M. Kubota, K. Ono, S. 
Komiya, Y. Ando, H. Eisaki, and S. Uchida: Phys. Rev. Lett. 

103 (2009) 037004. 

132) W. A. Groen, D. M. de Leeuw, and L. F. Feiner: Physica C 
165 (1990) 55. 

133) S. Ideta, K. Takashima, M. Hashimoto, T. Yoshida, A. Fu- 
jimori, H. Anzai, T. Fujita, Y. Nakashima, A. Ino, M. Arita, 

H. Namatamc, M. Taniguchi, K. Ono, M. Kubota, D. H. Lu, 
Z.-X. Shen, K. M. Kojima, and S. Uchida: Phys. Rev. Lett. 

104 (2010) 227001. 

134) K. Tanaka, W. S. Lee, D. H. Lu, A. Fujimori, T. Fujii, Risdiana, 

I. Terasaki, D. J. Scalapino, T. P. Devereaux, Z. Hussain, and 
Z.-X. Shen: Science 314 (2006) 1910. 

135) B. Lake, H. M. Ronnow, N. B. Christensen, G. Aeppli, K. 
Lcfmann, D. F. McMorrow, P. Vorderwisch, P. Smcibidl, N. 
Mangkorntong, T. Sasagawa, M. Nohara, H. Takagi, and T. 
E. Mason: Nature 415 (2002) 299. 

136) R. I. Miller, R. F. Kiefi, J. H. Brewer, J. E. Sonier, J. 
Chakhalian, S. Dunsiger, G. D. Morris, A. N. Price, D. A. Bonn, 
W. H. Hardy and R. Liang: Phys. Rev Lett. 88 (2002) 137002. 

137) J. M. Tranquada, B. J. Sternheb, J. D. Axe, Y. Nakamura, 
and S. Uchida: Nature 375 (1995) 561. 

138) Y. Sidis, C. Ulrich, P. Bourges, C. Bernhard, C. Niedermayer, 
L. P. Regnault, N. H. Andersen, and B. Keimer: Phys. Rev. 
Lett. 86 (2001) 4100. 



